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We consider a system of interacting fermions in two dimensions beyond the second-order per- 
turbation theory in the interaction. It is shown that the mass-shell singularities in the self-energy, 
arising already at the second order of the perturbation theory, manifest a non-perturbative effect: 
an interaction with the zero-sound mode. Resumming the perturbation theory for a weak, short- 
range interaction and accounting for a finite curvature of the fermion spectrum, we eliminate the 
singularities and obtain the results for the quasi-particle self-energy and the spectral function to all 
orders in the interaction with the zero-sound mode. A threshold for emission of zero-sound waves 
leads a non-monotonic variation of the self-energy with energy (or momentum) near the mass shell. 
Consequently, the spectral function has a kink-like feature. We also study in detail a non-analytic 
temperature dependence of the specific heat, C{T) oc T^ . It turns out that although the interaction 
with the collective mode results in an enhancement of the fermion self-energy, this interaction does 
not affect the non-analytic term in C(T) due to a subtle cancellation between the contributions from 
the real and imaginary parts of the self-energy. For a short-range and weak interaction, this implies 
that the second-order perturbation theory suffices to determine the non-analytic part of C{T). We 
also obtain a general form of the non-analytic term in C{T), valid for the case of a generic Fermi 
liquid, i.e., beyond the perturbation theory. 

PACS numbers: 71.10.Ay,71.10.Pm 



I. INTRODUCTION 

The validity of the Landau Fermi liquid (FL) theory continues to be a subject of intense discussions over the last 
five decades. In essence, the FL theory states that the behavior of interacting fermions at the lowest energies is 
similar to that of non- interacting fermions In particular, specific heat, C {T) , scales linearly with temperature 

T at T ^ 0, whereas the spin susceptibility, x (T) , approaches a finite value at T — > 0. This theory has been 
enormously successful in describing He"^ and a large number of metals. Yet, it fails to explain the properties of high- 
temperature superconductors Q and quite a few heavy- fermion materials . Given that the FL theory is essentially 
a phenomenological one, built on a number of appealing albeit unproven assumptions, the interest in its relation to 
microscopic models, which allow for controllable perturbative treatment, has always been intense since the times of the 
FL's inception. Perturbative calculations of the 1950s within a model of a low-density 3D Fermi gas with a repulsive, 
short-range interaction (a "non-ideal Fermi gas model" 0,0) reproduced the Fermi- liquid results for thermodynamic 
quantities up to leading terms in parameter T/Ep, where Ep is the Fermi energy (we set h ~ kB — I throughout 
the paper). These calculations were later extended to arbitrary dimensionality D with the result that the FL theory 
is valid for D > 1, provided that the interaction falls off with the distance rapidly enough, but becomes invalid in 
one dimension. The energy range for the FL theory, however, may shrink as the system approaches a quantum phase 
transition of some kind ^ . 

At the same time, perturbative calculations show that the similarity between a Fermi liquid and ideal Fermi gas does 
not go beyond the leading order in T / Ep. In a Fermi gas, sub-leading terms in C (T) /T and x (T) form regular series 
in powers of T^. For an interacting system, the second-order perturbation theory shows that already next-to- leading 
terms in T/Ep are non-analytic in D < 3. In Z? = 3, the sub-leading term in C (T) /T is T^lnT 6,, ,7, & 9, JLfl]. 
In addition, a non-uniform spin susceptibility, x((5), scales as Q^\nQ, where Q is the boson momentum [ll|. The 
non-analytic behavior becomes more pronounced as the system dimensionality is reduced. In 2D, the sub-leading 
terms in C (T) /T and in x (Q, T) scale as T [ll[ll[ll and max{Q, T} |Tl|ll|ll[ll[l3, respectively Non-analytic 



2 



corrections to the FL behavior have been observed in a number of experiments • Quite recently, the interest to 
these corrections has been revitahzed due to their importance for the effective theories of quantum critical phenomena 
in itinerant ferromagnets 0, 0, of the Hertz-Millis-Moriya type . 

The non-analytic behavior of the thermodynamic quantities is related to the long-range effective interaction between 
fermions which falls off as a power law of the inter-fermion separation, even if the nominal interaction is short-range. 
An example of such an interaction is scattering from the Friedel oscillation imposed by a static local perturbation, e.g., 
an impurity. This interaction falls off as with the distance from the impurity. Scattering from the r~^— potential 
of the Friedel oscillation leads to a non-analytic energy dependence of the scattering amplitude near the Fermi level 
and, as a result, to non-analytic corrections to the tunneling density of states 0| and conductivity ^ disorder- 

free-system, fermions interact via dynamic charge- and spin-density fluctuations. The non-analyticities in the dynamic 
density-density correlation function (polarization bubble) also give rise to a long-range retarded interaction falling off 
as r~^^~^^ These non-analyticities are due to the processes with both small and large (2fci?) momentum 

transfers, where kp is the Fermi momentum. 

The second order perturbative analysis in 2D not only reveals a non-analytic behavior of C(T) and x((5,T), but 
also brings about an unexpected result. Namely, the imaginary part of the fermion self-energy diverges logarithmically 
on the mass shell ui = ek |23, 24,.25j, J.3], if one linearizes the single-particle spectrum, e^, near the Fermi level. This 
log-singularity is the 2D analog of a stronger-power-law-singularity in ID ("infrared catastrophe") 26] . Although the 
mass-shell divergence does not affect the specific heat and the spin susceptibility to second order in the interaction, 
it does signal a potential breakdown of the perturbation theory in 2D. 

The purpose of this paper is to analyze the non-analytic corrections to the FL behavior beyond the second-order 
perturbation theory. Specifically, we focus on two issues. The first one is what happens to the mass-shell singularity 
beyond second order. Power counting shows that the mass-shell singularities proliferate with the order of the per- 
turbation theory. At first glance, this confirms Anderson's conjecture that the FL is destroyed in 2D j23|. However, 
this issue can be addressed properly only after a re-summation of the perturbation theory, which is what we will do 
here. The second issue is whether collective modes, which emerge once the perturbation theory is summed up to 
all orders, give rise to extra non-analytic corrections to thermodynamic quantities. The role of collective modes - 
zero sound for neutral fermions and plasmon for electrons - is especially intriguing for D = 2. In this case, power 
counting combined with the assumption that a collective mode is a free excitation (similar to a phonon) shows that 
the collective mode contribution to C (T) /T scales as T, i.e., it has the same form as the perturbative correction [T^ . 
This argument needs to be treated with caution, however. Indeed, since a collective mode arises at the infinite order 
in the interaction between fermions, it is unclear whether it can be treated as a free boson mode. We will see that the 
issues of the collective-mode contribution to C(T) and mass-shell singularities in the self-energy are related; namely, 
a contribution to the specific heat from the collective mode can be viewed as coming from the self-energy obtained 
by summing up mass-shell singularities to all orders in the perturbation theory. 

This paper is organized as follows. In Sec. |n]we introduce relevant scattering processes. In Sec. IIIII we discuss 
the self-energy of 2D fermions with both contact and finite-range interactions. In Sec. IIII Al and IIII Bl we analyze 
the mass-shell singularities in the imaginary part of the self-energy arising in the order-by-order perturbation theory. 
Re-summation of the perturbation series for the vertex part is performed in Sec. IIII CI The imaginary and real parts 
of the self-energy upon re-summation are discussed in Sees. 1111*01 and IIII El respectively. In Sec. IIII Fl we demonstrate 
that the spectral function exhibits a non-monotonic variation near the mass shell due to the interaction of fermions 
with the zero-sound mode. In Sec. IIVI we discuss the non-analytic contribution to the specific heat in two ways. 
First, in Sec. IIV Al we find C (T) via the self-energy, utilizing the results of Sec. IIIII Then, in Sec. II V Bl we evaluate 
the the non-analytic part of the thermodynamic potential directly. In Sec. IIV CI we consider the specific heat in a 
generic Fermi liquid. In Sec. IIV C 31 we consider the case of a Coulomb potential. Our conclusions are given in Sec.lVl 
Details of some of the calculations are presented in Appendices A-F. 

For the convenience of a reader, we present below a summary of the main results of this paper. 

A. Summary of the results 

1. Self-energy 

In Sees. IIIIIIII we consider mostly a 2D system of fermions with a weak, short-range repulsive interaction, specified 
by its Fourier-transform U{Q). The self-energy of such a system consists of two parts: a analytic one and a non-analytic 
one. The analytic part of the self-energy. 



San = A- be^ + c iuj'^ , 



(1.1) 
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where a, 6, and c are real, is determined by scattering events with large momentum transfers, of order kp. In this 
paper, we will be interested only in the non-analytic part of the self-energy, which comes from two types of effectively 
ID scattering processes. In the first type ("forward scattering"), all four momenta-two incoming and two outgoing- 
align almost along the same direction [cf. Fig. ^a)]. In the second type ( "backscattering" ) , both the initial and 
final momenta of the fcrmion pair are close to zero, while the momentum transfer can be near either zero or 2kp 
[cf. Figs.n^b) and (c)]. The angular spreading of the trajectories shrinks in the low-energy limit in proportion to 
|a;| /Ep for both types of scattering. To second order in the interaction and for a linearized single-particle spectrum 
[e/c = vp(k — kp)], the forward (F) and backscattering (B) contributions to the self-energy near the mass shell 
are [Ullllli, ^ 

Rel]«(c. = £,) = 0, Iml]«(^,fc) = ^ ^ In^^- (1.2a) 

ReE^iu. - = 8 IniSB(^,fc) = ^ ^ In^. (1.2b) 

Here Ep — kpvp/2 is the Fermi energy, 

mf/(0) mU{2kp) . ^, 

-0=^^, "2fc. = ^^ (1.3) 

are the dimensionless coupling constants which are assumed to be small, and 

A = cj-efc, (1.4) 

is the "distance" to the mass shell. On the Fermi surface (e^ — 0), ImS^ = ImEp + ImSg reduces to a familiar form 
ImE^(u;, kp) cx u? In |w| 

The special role of backscattering processes for the non-analytic corrections to thermodynamic variables of a FL 
has been considered earlier by two of us p^. In this paper, we present a complete description of forward-scattering 
processes. The peculiarities of these processes show up already at the second order: we see from Eq. 11.2all that on the 
mass shell, where A = 0, ImSp diverges logarithmically. The divergence is regularized 113,13,12^, |13| bv restoring a 
finite curvature of the single-particle spectrum, = d^ek/dk\, where kj_ is the component of k transverse to the 
local Fermi velocity Vi?(k). Finite curvature brings in a new scale 

Ac = Lj'^/W, (1.5) 

where W = mcVp/2, and the logarithmic singularity in Eq. H1.2a|) is rounded off at A ~ Ac. 

At a first glance, this regularization stabilizes the perturbation theory. However, starting from the third order in 
the interaction, the divergences due to forward scattering become of a power-law form; for a linearized spectrum, we 
find 

ImSf (u;,fc)oc (u^|c^/A|)"/2-i, (1.6) 

where n is the order of the perturbation theory. Finite curvature rounds off the power-law singularity on a scale A ~ Ac 
at every given order, but the resulting series for the self-energy holds in parameter Wc/I^j , where uJc = UqW/2, and does 
not converge for |cl;| < uJc- Therefore, the perturbation theory in 2D must be re-summed even for an infinitesimally weak 
interaction and realistic fermion spectrum. We show that the most divergent contributions to the forward-scattering 
part of the self-energy can be re-summed exactly to all orders in uq, without exploiting RPA-type approximations. 

Upon re-summation, the origin of the mass-shell singularities in the perturbation theory becomes clear: they 
correspond to the interaction between fermions and the zero-sound (ZS) collective mode. At every finite order of 
the perturbation theory, the collective mode coincides with the upper edge of the particle-hole continuum, and this 
degeneracy generates divergences in Ep . Once the perturbations are summed up to all orders, the ZS mode splits off 
from the continuum, and the power-law divergences disappear. The remaining logarithmic singularity is eliminated 
by the finite curvature. 

The total self-energy after re-summation is described by a scaling function of two variables A/ Ac and A/ A*, where 
Ac is defined in Eq. (|1.5() . and 

A* = uIlj/2. (1.7) 

is the scale at which perturbation series Eq. (|1.6|l diverges for a linearized spectrum. A general form of the scaling 
function is rather complicated, and will be discussed in the main text of the paper. In the limit A — > 0, when both 
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scaling variables are small, and also for low frequencies, reduces to 

Iml]«(. = (1.8b) 
Att hp \uj\ 

Comparing Eq. (|1.8a|) . (jl.Sbp with the second-order self-energy, [Eq. (|1.2a|l . Ijl.2bp ]. we see that the re-summation (i) 
eliminates the divergence in Iml]^(w = e^) and (ii) drastically changes the result for KeY^^{uj = e^). In particular, for 
a constant interaction {uq = U2kF ) , re-summation of higher-order terms in the self-energy cancels out the second-order 
term, so that full ReS^(ci; — ek) vanishes. 

Perhaps the most essential result of our analysis is that the interaction with the zero-sound not only leads to 
a scaling behavior of the self-energy, but also results in a singularity of the self-energy: the derivative dhnYi^ / d/S. 
diverges as l/VA — A* at A = A*. This singularity is present in the non-perturbative regime, i.e., for \bj\ < ujc. 
Physically, it corresponds to a change in kinematics of ZS waves emission. On the mass-shell (A = 0), emission of ZS 
waves by fermions is impossible as the zero-sound velocity is larger than the Fermi one. For < A < A*, emission is 
possible but it is subject to a Cherenkov-type restriction: a fermion with frequency lu > can only emit a ZS wave 
in the frequency interval fl < u>A/A*. For A > A*, emission of ZS waves in the whole interval < f2 < w becomes 
possible. The self-energy is singular right at the onset of the Cherenkov-type restriction for emission of ZS waves. 

The singularity in the self-energy translates into a kink in the spectral function, A{uj,k) — —7T^^lmG^{uj, k) at 
A = A* (cf. Fig. IS]). This effect is actually present for both short-range and Coulomb interaction. For the latter, 
the collective mode is a plasmon, and the kink is positioned near the Fermi surface, where ek — (or A — lu). 
The prediction for a kink A(uj, k) can be verified in angle-resolved photoemission on layered compounds 29] or in 
momentum-conserving tunneling between two parallel layers of 2D electron gas ,30] . 



2. specific heat 



In the second part of the paper fSec. lIVI) . we analyze the non-analytic behavior of the specific heat, C(T), for three 
types of interaction: i) short-range, weak repulsion; ii) Coulomb interaction; and iii) generic Fermi-liquid interaction. 
For all three cases, we find that a non-analytic term in C(T) behaves as T^. In the perturbation theory, an origin of 
this term can be simply related to a non-analytic, form of the real part of the self-energy. To second order in 
the interaction, a non-analytic part of C(T) 



SCiT) = C{T) - 7T, 

where 7 is the (interaction-dependent) Sommerfeld factor, was shown earlier |l3l | to be 

9C(3), 



5C (T) /T : 



UoU2kf 



-Cfg/Ef, 



(1.9) 



(1.10) 



where 



Cfg 



TO7rT/3 (1.11) 

is the specific heat of a Fermi gas [Eq. (|l.ll|l ] (see also Refs. '3Q|^). 

The issue considered in this paper is whether SC{T) is affected by the interactions of fermions with the zero-sound 
mode. At a first glance, it should be. Indeed, C(T) is related to an exact retarded Green's fimction {ut, k) — 

[t^-efe + S«(t^,fc)]"^ via Q 



2 d 

C[T)IT = — — 

TT 01 



dPk 

(2^ 



dujuj^—^ arg G^'(u!, k) 

OUJ 



(1.12) 



where np is the Fermi distribution function. As the real part of the self-energy is changed significantly by a non- 
perturbative contribution from the ZS mode, the corresponding change in G^ should a priori affect G{T). However, 
we see from Eq. (I1.12f) that-contrary to the common wisdom- not only the real but also the imaginary part of T,^ 
affect C{T). Indeed, at low frequencies, which we only need at small T, both perturbative and non-perturbative parts 
of the self-energy are asymptotically smaller than \u}\; thus in Eq. (|1.12|l can be expanded to first order in 
with the result 



2 d 




r (Pk 




n df 


T J 


1 (27r)2 


J — 00 



dujLu^ {ReS]^(w, fc)ImG^(w, k) + ImE^(w, fc)ReG^(w, fc)} 

ULU 



(1.13) 
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Substituting the second-order result for ReE-'^, Eq. (|1.2b|) . into Eq. (|1.13|l . we indeed reproduce the T-dependence of 
SC{T)/T, as given by Eq. (|1.1U|) (more care is required to reproduce a numerical prefactor as it turns out that one 
should use an expression for ReT,^{uj, k) at finite temperatures- cf. Sec. lIV Ajl . The perturbative part of lui>l^{uo, k) 
does not contribute to the specific heat as it depends on uj but not on fc; as a result, the second term in Eq. van- 
ishes by parity upon switching from integration over cPk to that over rfefe. However, the non-perturbative contribution 
to ImE^(Li;, k) due to the interaction with the zero-sound mode depends strongly on fc. As a result, both ReS-'^(cL', fc) 
and ImS^(a;,fc) contribute to C (T) . We show that non-perturbative terms in these two contributions cancel each 
other, i.e., there is no non-perturbative contribution to the specific heat. A non-analytic, T^-term in 6C{T) then 
comes entirely from the perturbative part of E^, and Eq. (|f .10() is the complete result for dC{T) to second order in 
the interaction. 

Another way to understand an absence of the non-perturbative contribution to the specific heat is to evaluate the 
thermodynamic potential, S, directly in Matsubara frequencies, and then use the relation between S and C (T) . This 
is done in Sec. IIVBI In contrast to the real-frequency description of the self-energy, there are no singularities at 
any order of the perturbation theory for S in Matsubara frequencies. This means that, as long as the interaction is 
weak, one can truncate the perturbative series at an arbitrary order and be sure that the higher-order terms would 
give only sub-leading contributions. In this approach, there simply cannot be non-perturbative contributions to the 
thermodynamic potential, and hence to C (T). To second order, this procedure gives the same result as in Eq. H1.10|l . 
For completeness, we also evaluate thermodynamic potential in real frequencies and demonstrate explicitly how the 
collective-mode contribution to the specific heat cancels out. 

Finally, we extend our analysis to a Fermi liquid with not necessarily weak interaction. We find that the term 
in the specific heat for a generic Fermi liquid is expressed via the charge (c) and spin (s) components, fc,s{Q)i of 
the quasi-particle scattering amplitude between particles at the Fermi surface, at angle 6 — tt between two incoming 
momenta as 

SCiT)/T = --^7^ W + 3/,^(7r)] T, (1.14) 

where Vp = kp/m* and m* is the renormalized effective mass. We remind the reader that the scatteringamplitude 
(as a tensor in the spin space) is related to a particular limiting form of the interaction vertex, r^(0), as [l|,l2 

f (6) ^ Z'r'^ (9) , (1.15) 

where 

t''{9)= lim f(fcFni,0;fcFri2,0|Q,f7), (1.16) 

where r(ki, o^i; k2, W2IQ, ^) is the vertex for a process (ki, wi; k2, 0^2 ^ ki — Q, a; — £7; k2 Q, -I- fi), and 6 is the 
angle between ki and k2 . In the FL theory, the renormalizations of the thermodynamic quantities are expressed via 
the angular moments of the Landau interaction function, which is related to another limiting form of the vertex, 

f'^(e') = lim f(fci.fii,0;fcj.n2,0|Q,17). (1.17) 
Q/|r2| *o 

Simple algebraic relations between the partial components of F*^ (9) and F"^ {9) enable one to express the analytic parts 
of thermodynamic quantities either via the moments of the Landau interaction function or that of the scattering 
amplitude. However, the non-analytic, T^- part of the specific heat is related to r'^(0) at a particular angle {9 = tt), 
rather than to its angular average. As there is no simple relation between T^{9) and T'^{9) for any given angle, 
including 9 — 7:, the non-analytic in the specific heat in general cannot be expressed in a compact form in terms 
of the Landau interaction function without making additional model approximations @. In this respect, our result 
for SC{T) differs from that of Ref. 31], where 6C{T) was expressed via the charge and spin components of a single 
Landau parameter F'~' (x J d9T'^{9). We did, however, obtain 6C{T) in terms of r'^(7r) in the limit when its charge 
component is much larger than the spin one. The limit when the charge component tends to infinity whereas the spin 
one tends to zero describes the Coulomb interaction in the high-density limit. In this case, we find that the term 
in C(r) is universal and independent of the electron charge [cf. Eq. ()4.57|l ]. in agreement with Ref. jsj- 
In the rest of the paper we present the details of our analysis. 



II. SCATTERING PROCESSES 



In a typical event of interaction between low-energy quasi-particles with momenta fci w fc2 ~ kp, the change in 
the momentum of a given quasi-particle, Sk = |ki — k'^^l, is of order kp, but not necessarily close either to zero or to 
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FIG. 1: Scattering processes responsible for divergent and/or non-analytic corrections to the self-energy in 2D. a) "Forward 
scattering" -an analog of the "(J4" -process in ID. All four fermion momenta are close to each other, b) Backscattering-an analog 
of the "(;2"-process in ID. The net momentum before and after collision is small. Initial momenta are close to the hnal ones. 
Although the momentum transfer in such a process is small, we still refer to this process as "backscattering" (see the discussion 
in the main text), c) Another component of the backscattering process: 2kF— scattering. 

2kp. These large-angle scattering events are responsible for the analytic part of the self-energy, Eq. Hl.l|l . In addition, 
there are special scattering events in which either Sk ~ |ri|/wF ^ kpOT \5k — 2fci?| ~ |51|/w_f kp, where il is the 
energy transfer. Although the phase space associated with these events is small for D > 1, these processes give rise 
to non-analyticities in the dynamic density-density correlation function and eventually determine non-analyticities 
in S(w) [iJI- The role of these special processes increases as the dimensionality is reduced. For D = 3, the non- 
analytic part of the self-energy E,ia oc uj^ln(—iuj), resulting from the special processes, is sub-leading to the analytic 
one, resulting from the generic processes. However, already for D = 2, the non-analytic part (S)^^ oc ioj'^ In (— iw) 
dominates over the analytic one. 

In 2D, kinematics of processes with small momentum is essentially one- dimensional, i.e., the initial and final 
momenta of two interacting fermions are either almost parallel or antiparallel to each other. (In 3D, both ID and 
non-ID processes contribute to the non-analytic behavior.) Accordingly, these processes can be divided into two 
types. In the first type, the two colliding particles move initially almost in the same direction (ki « k2) and retain 
their respective momenta after the collision, so that all four momenta (two initial and two final) are close to each 
other 

k'l « ki « k2 « ka'. (2.1) 

This type of process is shown in Fig. ^a). In g-ology [s^, such an event is called "g4-scattering" . The deviation 
from the purely ID kinematics is due to finite energy transfers: a typical angle between momenta in Fig. ^a) is of 
order \uj\/vpkp <C 1. In what follows, we will refer to the process in Fig.^a) as simply "forward scattering". In the 
second type, the colliding particles move initially in almost opposite directions (ki sa — k2) but, as for the forward 
scattering case, they also retain their respective momenta after the collision. The difference between such an event 
and the forward-scattering one is that not only the transferred but also the total initial and final momenta are small. 
This type of process is depicted in Fig.^b). In 5— ology notations, this is a "g2-process" . 

Another process which contributes to the non-analytic part of S^(cj) is a "2kp- process" (or f;i-scattering, in 
f/-ology notations), in which two fermions moving initially in almost opposite directions, reverse their respective 
momenta [see Fig. ^c)]. Since both processes in Figs, ^b) and (c) contribute to the same scattering amplitude 
/(ki,k2) with the angle between initial momenta ki and k2 being close to tt, we will refer to both of them as 
"backscattering". To distinguish between the two, we will refer to Fig. ^b) as "(72— backscattering" and to Fig.^c) 
as "2kp— backscattering". 

Scattering by 2kp is one-dimensional in all dimensions. As we just said, forward- and backscattering become one- 
dimensional in _D = 2. Wc thus conclude that for D = 2 the non-analytic part of the self-energy comes from essentially 
ID scattering processes, embedded into the 2D phase space. 

We pause here for an important remark. Although there is a strong similarity between special scattering processes, 
resulting in non-analytic behavior in 2D, and their ID analogs, there is also an important difference. Namely, neither 
f/2 nor 54 processes lead to a non-analytic behavior of thermodynamic quantities in ID. This is already obvious from 
the fact that a ID Hamiltonian with a linearized spectrum and in the absence of gi-scattering (Tomonaga-Luttinger 
model) allows for an exact diagonalization in terms of new excitations-free bosons. As a result, the specific heat is 
strictly linear in T and the spin susceptibility is simply a constant within the Tomonaga-Luttinger model. However, 
if gi-scattering is present even as a marginally relevant perturbation, exact diagonalization in terms of free bosons is 
no longer possible, as the spin sector is now described by the sine-Gordon rather than Gaussian theory. This results 
in strong non-analyticities in both the specific heat dC(T) (oc TlnT) and spin susceptibility dxs (oc | IniJj), where H 
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a) 




FIG. 2: Non-trivial second (a) and third (b) order diagrams for the self-energy. 

is the magnetic field It is possible to obtain these ID non-analyticities within the same frajmework as their 

2D analogs are analyzed in this paper, but we defer this discussion to a separate publication |35ij. 

III. SELF-ENERGY 

In this Section, we derive an expression for the self-energy to all orders in the interaction. We analyze and re- 
sum the mass-shell singularities in the forward-scattering part of the self-energy, and also review the behavior of the 
backscattering part. For the sake of completeness, however, we start with the brief discussion of the second- order 
results for the self-energy. The fermion self-energy is defined via the Dyson equation 

G-i=G(7i+S, (3.1) 

where Gq and G are the bare and exact Green's functions, respectively. (Notice that we define E with an opposite 
sign compared to Refs. 

A. second order 

To second order in the interaction, there are only two non-trivial diagrams for S, which are shown in Fig.|2Ia). For 
a contact interaction, V (r) = US{r) , the contribution from diagram (a2) is (—1/2) of that from diagram (al). The 
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net contribution to the imaginary part of the retarded self-energy at T = is given by 

— / — ^ImG^ (w + r!, k + Q) Imn^(r!, Q), (3.2) 

-u J (27r) 

where (w, k) = [ui — ek + jO+) ^ is the free retarded Green's function, n^(ri, Q) is the polarization bubble of free 
fcrmions, and subindex 2 of the self-energy denotes the order of the perturbation theory. A general expression for 
Il^{rL^ Q) is rather complicated 36], but in what follows we will need only its two asymptotic forms. The first of these 
forms is valid for small Q 



n-(f7, - -i 1 1 + / I ' (3-3) 



and the other one is valid near Q = 2kp 



m 




Tl^i^,Q) = -^\^-\^^^ + \l[^^^] -['-i^T—] I I- (3.4) 

Here m is the fermion's mass and is the Fermi velocity. Non-analyticities in the two limiting forms of the bubble 
describe Landau damping and Kohn anomaly, respectively. Landau damping of an excitation with energy fl and 
momentum Q is possible only within the particle-hole (PH) continuum, i.e., for CI < vpQ. For CI <^ vpQ, the non- 
analytic term in Eq. 1)3. 3|l scales as n^jjg(r2, Q) oc iCl/\Q\. The Kohn anomaly near 2kF is static for Q > 2kp [in this 
range, the non-analytic part of behaves as n^jjg(0, Q) oc {Q — 2kpy^'^] but is dynamic for Q < 2kp [in this range, 
n,tg(f7,Q)cxzr!/(2fcj.-Q)i/2]. 



1. hacks cattering 

To logarithmic accuracy, both 52— and 2fci^— processes [Fig. ^b) and (c), respectively] contribute equally to the 
non-analytic part of the fermion self-energy. For a contact interaction, the sum of the two contributions is |13j 

ImS«B(-,fc) = ^|^ In^, (3.5) 



where B stands for backscattering and u is defined in Eq. H1.3|l . It is important in what follows that, to logarithmic 
accuracy, ImE^g depends only on lo but not on efe. For the sake of completeness, we present the derivation of Eq. I|3.5|) 
in Appendix 1X1 We find that for 2A;j;^— processes, a non-analytic part of the self-energy originates from the dynamic 
Kohn anomaly, whereas the static Kohn anomaly contributes only to the regular part. 

Higher-order contributions to the backscattering part of the self-energy form regular series in u which result in 
the renormalization of the prefactor. As we keep u small, it suffices to stop the perturbation theory at order u^. 
Consequently, we set 

Eg = E«B (3.6) 

in the rest of the paper. 



2. forward scattering 

The second-order forward-scattering contribution to the self-energy is given by |13| 

where F stands for forward scattering, and we remind that A = w — e/c is the "distance" to the mass shell. Away from 
the mass shell, this contribution behaves as w^lnjw], i.e., it has the same functional form as ImE5. In contrast to 
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the backscattering part, however, the forward-scattering contribution diverges at the mass shell, i.e., for A — > 0. The 
origin of this divergence can be traced back to the form of the polarization bubble at small momenta. From Eq. 
we find that 

Imn^(»,Q)^-(^) ^L^ 9{vFQ~\n\), (3.8) 

where 9 (x) is the step function. Imll^ {rt,Q) has square-root singularities at = vpQ. On the other hand, 
expanding ek-i-Q in G^{uo + 17, k -I- Q) in Eq. H3.2|l as ek+Q = Cfc + vfQ cos 9 and integrating over 0, we obtain another 
square-root singularity 



/ 



d9lmG^ = (3.9) 



On the mass shell (w = e^), the arguments of the square roots in Eq. H3.8|) and Eq. (|3.9|) coincide, hence the integral 
over cPQ diverges logarithmically. DM H3.7(l is valid only for a linearized fermion dispersion. Two of us demonstrated 
in Ref. |l3l | that finite curvature of the dispersion eliminates the logarithmic singularity in ImS^p. To keep our 
presentation uninterrupted, we continue to proceed with the analysis of the singularities due to forward scattering, 
assuming that the curvature is equal to zero, i.e., the dispersion is linear. We then discuss separately the modifications 
imposed by a finite curvature of the dispersion (cf. Sec. lIIID4)l . We emphasize again that the logarithmic singularity 
in the self-energy arises from essentially ID scattering processes, embedded in a 2D phase space. Therefore, this 
sing ularity can be viewed as a pre-cursor of a stronger (power-law) singularity in ID ("infrared catastrophe") [2^. 
Ilia. 



B. higher-order forward-scattering contributions 



Higher orders of the perturbation theory contain more bubbles with small momenta ("soft bubbles"). As a result, 
the mass-shell singularities proliferate. The third-order diagrams, shown in Fig. |21 contain the square of the soft 
bubbles. These bubbles appear either explicitly (as in diagram bl) or are generated upon integrating over fermion 
energies /moment a in the rest of the diagrams. The singular part of the self-energy at this order is given by 

^ / -^ImGo«(c. + f7,k+Q)n^(f7,g), 
-UJ J (27r) 



where 



(17, Q) EE Imtf + iO+, Q) . 
The most singular term in 11^(17, Q) is given by 

[n^(r!, Q)] - - {m/2^f Tin \n\ 6 {n^ ~ vIq^) 



(3.10) 



(3.11) 



The product of the square-root and delta-function singularities (from J G'^c?6' and Il§{n,Q), respectively) gives rise 
to a one-sided, square-root singularity on the mass shell: 



ImEf F (^^, fc) = - 



\/2u3 



(I) 



20 EfM A 
It can be readily verified that at n— th order 

for n > 2. Collecting forward-scattering contributions to all orders in w, we obtain 



Ef ^ / 2 ^ ^ 



n=l 



(3.12) 



(3.13) 



(3.14) 



where C„ are the numerical coefficients. We see that perturbative expansion in u works only for u^uo/ IS. <^ 1. Outside 
this range, series in u does not converge, and one needs to re-sum the perturbation theory. 
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C. Re-summation of forward-scattering contributions 



To perform the re-summation of the perturbation theory, we need to select diagrams with the maximum number 
of particle-hole bubbles at small frequency/momentum. It is convenient to select first analogous diagrams for the 
four-fermion vertex, Tq,^ ,yg(pij?;p3,p4), and then relate E to F via the Dyson equation. In this subsection, we will be 
using notations p = k) and q = {Q.m,Q)-, where = 7r(2m -|- 1)T and flm = 2i:mT. 



1. four-fermion vertex, zero-sound mode 



The diagrams for T with the maximum number of particle- hole bubbles form familiar ladder series (see Fig. IJJ, 
when F is anti-symmetrized with respect to a permutation of either initial or final states. However, the procedure of 
finding an overall prefactor at order v is somewhat involved 0, as it requires counting the number of diagrams at 
the same order in a conventional diagrammatic technique operating with a non-symmetrized vertex, F. We choose to 
sum the diagrams for a non-symmetrized vertex to all orders first, and then anti-symmetrize the result. The second- 
and third-order diagrams for F are shown in Fig. |3 A general procedure of summing such diagrams to all orders is 
described in Appendix El K leads to the following result for f 



ra/3,7e(pi,P2;Pl " 9,^2 + = r((?) = -U 



where ct°o are Pauli matrices (a = x,y,z), and 



1 



c/n {q) ■■ 



1 



2 1 + [/n (q) 



(3.15) 



(3.16) 



are the (dimensionless) charge- and spin vertices, respectively. An anti-symmetrized vertex is obtained from F by the 
following procedure 



(3.17) 



For the case of t/ > 0, which we are interested in, the retarded charge vertex, ^ has a pole determined from the 
equation 1 — UTi^{ci, fi) = 0. A two-particle excitation corresponding to the pole in Qp is a zero-sound collective 
mode. Since 11^ is real for fi^ > VpQ'^, and can be arbitrarily large (and positive) when \n\ approaches vpQ, the 
zero-sound pole exists already for an arbitrarily small U. Near the pole, Gp is of the form 



Qp 



(f7 + iO+) - c2Q2 



where c is the zero-sound velocity 



c = vfVI + (1 + 2u) ~ Wi^(l + M^/2) > vp. 



(3.18) 



(3.19) 



We see that zero-sound velocity c differs from vp only by a w^-term. This means that the zero-sound mode = cQ 
is just above the upper boundary of the particle-hole continuum, which, for small Q, is given by vfQ- We also see 
from Eq. H3.18|l that the quanta of zero sound are not free bosons as the residue of the zero-sound pole in Eq. (|3.18|l 
is proportional to . The spin vertex {Q^) also has a pole, but it is located on the imaginary axis. Consequently, the 
corresponding collective mode is over-damped. As an independent check, we verified that the diagrams in Fig. |3 sum 
up to zero for the case of spinless fermions. This result is a manifestation of the Pauli principle: spinless fermions do 
not interact via contact forces as the Pauli principle forbids them to be at the same point in space. 



2. Dyson equation 

The self-energy due to forward scattering is related to the vertex function via the Dyson equation Q 

^F,M = f UG{p-q) ~ f UT^^,^p{p,p' +p" -p-p\p")G{p')G{p")G{p' +p" ^p), (3.20) 
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a 



► s ► 




Y 



► ^ ► 






FIG. 3: Vertex diagrams with maximum number of particle-hole bubbles to third in the interaction, 
obtained from those in the second column by a permutation a — > /3, 7 — > e, are not shown. 



Additional diagrams, 



where 



(3.21) 



In principle, the Green's functions in the Dyson equation are the exact ones. However, it can be verified that self- 
energy insertions into the diagrams diverging near the mass shell do not give rise to additional mass-shell singularities. 
As we keep u small, regular corrections are thus irrelevant, and we can safely use bare G"s instead of the exact ones 
in Eq. (|3.20(l . Substituting Eq. H3.15|l into Eq. H3.20|) . we obtain T,F,ai3 = Sap'Sp where 



U + U'Uiq) + - 



1 U'^lP{q) 3 U^n^iq) 



2{i-un{q)) 2(i + [/n(g))J 



G{p-q). 



(3.22) 



Technical details of the derivation leading to Eq. H3.22|l are presented in Appendix ^ Expanding Eq. H3.22|l to third 
order in U, we reproduce the results of the conventional perturbation theory, Eqs. (|3.7|l and (|3.12|l . We remind that 
in the perturbation theory ImS^ diverges upon approaching the mass shell- logarithmically to second order in U, 
and as l/-^a; — Ck to third order. It is convenient to rearrange the terms in (|3.22|) and decompose Sp into three parts 
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making use of the charge- and spin vertices, introduced in Eq. (|3.16() . as 

Sf(p) = Sp (p) + (p) + Ecx; (3.23a) 
Sp(p) = U I gp{q)G{p~q)- (3.23b) 



E,(p) = W I g„{q)G{p-q); (3.23c) 

Sex = / [2C/ - f/2n(g)] G (p - g) . (3.23d) 

Terms Sp and So- correspond to the interaction in the charge- and spin channels, respectively, and are summed to all 
orders in U. The remainder, Sgxi contains extra contributions of the first and second orders in ?7, not included in the 
first two terms, reproduces the second-order result Eq. 1)3. 7f) . 

Before wc proceed further, a comment is in order. Our results for the vertex and self-energy formally coincide with 
those found in the paramagnon (spin- fluctuation) model jSj (except for the remainder term, Eox, in Eq- l|3.23a|l which 
was neglected in Ref. [^). However, our results have been obtained in a more general approach. In the paramagnon 
model, the self-energy is given only by diagrams of the type (bl) and (b5) in Fig. 13 i.e., it involves only RPA diagram 
in the charge-channel and ladder diagrams in the spin channel. We included all diagrams with the maximum number 
of bubbles and found that the overall combinatorial coefficients at each order are such that the summation to all 
orders results in two independent geometric series-one for the charge channel, and the other for the spin channel. 
It does not mean, however, that we have obtained exact results for V and S. Indeed, we considered only forward 
scattering, kept u small, and neglected all diagrams that constitute regular series in u and does not give rise to 
proliferating mass-shell singularities in S. From this perspective, the controlling parameter for our approximation is 
not the coupling constant u itself but a combined parameter u'^|aj|/|A| which measures the proximity to the mass 
shell. We sum up the series in u^|a;|/|A|, and neglect regular corrections in u at every order . 

D. Imaginary part of the self-energy to all orders in the interaction 

We now evaluate the forward-scattering part of ImS^ near the mass shell. The imaginary part of the retarded 
self-energy comes from two sources: from the particle-hole continuum (|J7| < vpQ) , where Imll^ ^ 0, and from the 
collective mode at = cQ, where Gj^ has a pole. The spin-channel part of the self-energy, ImE^, comes only from 
the continuum, whereas the charge-channel part contains contributions from both the continuum and collective mode. 
Accordingly, the imaginary part of the total self-energy can be represented as 

ImEf = ImE^H + Im^zs + ImSfx> (3-24) 

where(. . . )ph and (. . . )zs stand for the particle-hole and zero-sound contributions, respectively: 

ImE^H = (imS^ + ImS^Opjj , ImEfg = (imS^O^g • (3.25) 

1. remainder term T^ex 

Comparing Eq. (|3.23d|l and Eq. we see that the remainder term in decomposition Eq. H3.23a|l . ImE^, is op- 

posite in sign and equal in magnitude to the second-order forward-scattering contribution to the self-energy, Eq. (|3.7|l : 

2. particle-hole contribution 

Term ImEpjj contains the imaginary parts of the retarded vertices in the charge- and spin channels: 

„ 1 C/Imn^ 

Im^f = - 5 ^. (3.27) 

2 (1 ^ [/Ren«)' + (C/Imn«)' 
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The first term in the denominator of Eq. (|3.27() can be replaced by unity because for q < 2kp^ Rell-'^ = — TO/27r, i.e. 
-UReU^ = M < 1. Substituting then Eq. (|T77jl into Eq. H3.23b|) and Eq. I|3.23c|) . we obtain 



rImG^(a) + f7,k + Q) 



Imll^ 



,2 ■ 



(3.28) 



Substituting Imll-" from Eq. and ImG^ from Eq. into Eq. 

contribution, we obtain after some algebra 



47r Ef 



In 



Gi 



2A 

u'^\io\ 



1 + (c/imn«)' 

and keeping only the forward scattering 



(3.29) 



where 



Gi{x) = 21n2-l/2 + ln| 



2 Re 



zdz 



1 + ^1 - x/z 

yjl — xj Z 1 — ^1 — xj Z 



In 



(3.30) 



Subscript / in G/ implies that this is a scaling function for the imaginary part of the self-energy. We subtracted off a 
constant term in G/ so that G/(0) = 0. This is equivalent to neglecting a regular, cj^-contribution to ImEpjj. A plot 
of Gi(x) is presented in Fig. 01 At large and positive x, the integral term in scaling function Gi(x) falls of as a;"^/^, 
whereas for large and negative x, it falls off as (— x)~ . In either of these limits, G/ « In |a;|~^ and, consequently. 



47r Ef A 



Further expansion in powers of 1/x yields 



ImSl^pH 



,V2u3 ^2 

' 20 Ef 



(3.31) 



(3.32) 



In the opposite limit of small x, function Gi{x) vanishes as a; In |a;|. As a result, net ImSpjj remains finite at A = 0, 
and for a; <C 1 (i.e., for A <C u^uj) it behaves as 



An Ef 



In 



Ef 



2A ^ |A| 
■ in ■ 



(3.33) 



Comparing the limiting forms of Eq. I|3.31|l and Eq. (|3.33|l . we see that higher order terms in u simply cut the 
logarithmic divergence in ImSpjj for |A| < u^|cl;|. To logarithmic accuracy, one can then approximate ImSpjj by 



n v?' uj"^ Ef 

ImSpTT = -— m -; — r, 



(3.34) 



where w = max { A, it^|u;| }. The appearance of u under the logarithm in Eq. H3.34|l is a reminder that Eq. H3.34f) 
includes all orders of the perturbation theory. Indeed, a crossover between Eqs. (|3.31ll and (|3.33|l occurs at a scale 
|A| ~ it^|w|, which is not accessible within the perturbation theory. 



3. zero-sound contribution 



For the collective-mode contribution to the self-energy in the vicinity of the mass shell, i.e., for |A| <C |a;|, we obtain 
from Eq. H3.23b|l and Eq. (jXTfjl 



ImS" = 



4Tr 



-Re 



QdQ 



y^{2A/vFQ) sgntj - ^ 



An Ef' 



-Ft 



2A 



(3.35) 



For X > 0, scaling function Fj (x) is given by 



Fj (x) = 27r 



i{l,v^} 



dy 



371^x78, for X < 1; 
3|i,i,-i_^_V:F^(|x + l) 



for X > 1. 



(3.36) 
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FIG. 4: Scaling functions Gi (x) [ Eq. llOTl ] and Fi (x) [ Eq. iTM ]. 



For negative x, Fj(x) = 0. At x = 1, Fj{x) is continuous but its derivative is singular: dFj{x)/dx cx 1/y/x — 1 for 
a; > 1. A plot of is shown in Fig. ^ For x ^ I, i.e., for u ^ 0, 

« 2n/5y/x. (3.37) 

In this limit, Fj{x) cx u and thus ImS^g oc m'^, as is to be expected, as the zero-sound propagator, UGp, is of third 
order in U near the pole. Substituting this limiting form into Eq. II3.35II . we obtain for w/A > 

ImE^o = J — . (3.38) 

Combining Eq. (|3.38|) and the third-order particle- hole contribution as given in Eq. H3.32|l . we reproduce the result 
of the third-order perturbation theory, Eq. (|3.12() . Expanding (ImSp)^in powers of u further, we indeed reproduce 
the structure of higher order terms in the perturbation theory, Eq. I|j.l3|l . All these terms diverge when A — s- 0. 
However, the full result shows that the perturbative expansion in u for the zero-sound contribution is valid only for 
|A| ^ u^IljI. At |A| = u^|Li;|/2, i.e., at x = 1, ImS^g has a maximum. Upon further approach to the mass-shell, 
ImE^g decreases as A^ and eventually vanishes on the mass shell (A = 0). Vanishing of ImEfg on the mass shell 
is due to a Cherenkov-type restriction: because the zero-sound velocity c > vp, an on-shell fermion cannot emit a 
zero-sound boson, hence the fermion's lifetime becomes infinite. 

Observe that ImE^g is asymmetric with respect to a change in sign of A: ImS^g 7^ only if A and lo are of the 
same sign. This asymmetry follows simply from the energy and momentum conservation. For example, a fermion of 
energy w > above the Fermi level can emit a soft zero-sound boson of frequency < < a; provided that 

cj — efe = — vpQcosO = n[l — (vp/c) cos 6*] , (3.39) 

which, for c « (1 + is equivalent to 

< < A <2n <2uj. (3.40) 

Thus, emission is possible only if both uj and A are positive. A similar consideration for the case when a fermion of 
energy to < below the Fermi level absorbs a zero-sound boson of frequency fl in the interval (— 1^1,0) shows that 
absorption is possible only for A < 0. Eq. (|3.40l) also clarifies the meaning of a characteristic scale A* = uiu'^ /2. For 
A > A* emission of bosons with any frequency in the interval < J7 < w is possible. In particular, a fermion can 
emit only one boson of frequency Q = uj and "land" on the Fermi level. For A < A*, i.e., when a fermion is close 
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to the Fermi level, emission of bosons with frequency fl > 2A/m^ is impossible, and the fermion relaxes to the Fermi 
level via emitting a large number of low- frequency bosons. As a result, the relaxation slows down which corresponds 
to a decrease in ImEfs for A < A* . 

Combining the results for S^, Spjj and Sfs' ^'^^ that the summation of the power-law divergent diagrams for 
the self-energy leads to a non-trivial result: the total self-energy due to forward scattering undergoes a non-monotonic 
variation near the mass-shell. All power-law divergences of the form u"A^~"/^ are now eliminated. However, we still 
have a logarithmically divergent term ImS^, given by Eq. H3.26|l . As this term does not contain higher than the 
second order in u, its divergence can be cut only by a finite curvature of the fermion dispersion (see Sec. IIII D 4|) . 

Notice also that the collective-mode contribution to ImEp is smaller than the rest of the contributions by a large 
logarithm. Indeed, the maximum value of Fj{x) in Eq. (|3.35|) is of order one, so that ImS^g S u^i^j'^/Ep, whereas 

ImSf ~ ImS^jj ~ ImEf^ ~ (u^uj^/Ef) \tlEf/ > u^uo'^/Ef- (3.41) 

Still, ImS^g exhibits a non-monotonic and rapid variation near the mass shell at A ~ v^u) <C w, whereas other 
contributions are either constant or vary only smoothly on this scale. This feature will have consequences for the 
spectral function, discussed in Sec. IIII Fl 



4- modifications due to a finite curvature of the fermion dispersion 

As we have already mentioned in Sec. IIII Al the logarithmic mass-shell singularity in ImSp at the second order 
can be eliminated by accounting for the finite curvature of the dispersion. Technically, this amounts to retaining 
the quadratic-in-g term in the expansion of Ck-t-Q in Q. A straightforward analysis shows [l^ that the logarithmic 
singularity in the second-order diagram is cut at a certain distance to the mass shell |A| ~ Ac, where, we remind, 
Ac = u}'^ /W , W — kp/{2mc), and l/nic is the curvature. For the quadratic dispersion efe = fc^/(2m), rric = m, hence 
W = Ef- For a non-quadratic dispersion, W and Ef are not equivalent, but, at least for any power-law spectrum, 
they are of the same order. Therefore we will not distinguish between W and Ef in the rest of the paper. Cutting 
the log-singularity in ImS^p at Ac, we obtain 



ImI]^p(A = 0,c.) = f ^In^. 

An Ef \(^\ 



(3.42) 



The net second-order self-energy, i.e., the sum of backscattering and forward-scattering contributions, is then given 

by 



ImEf {uj) = ImEf [lo) -f ImEf p (w, A = 0) = 2ImEf {uj) ^ In 

27r Ef |w 



(3.43) 



The elimination of In |A| singularity at the second order does not eliminate the need for re-summation of the pertur- 
bation theory, since higher-order terms diverge as powers of |A|~^. Indeed, cutting the singularities at |A| = Ac in 
the general expression Eq. H3.14|) for ImSp (w), we obtain 



ImE^(w) 



9 "? 
U W 

47r Ef 



In 



Ef 



i/2 



where 



— „,2 



u^Ef 



(3.44) 



(3.45) 



Obviously, the series for ImEp does not converge for \bj\ < 



one still needs to re-sum the perturbation 



theory. We already know, however, that all power-law divergences are eliminated after such re-summation even for 
infinite ujc- Finite curvature is not going to modify the results for A 3> Ac. For arbitrary A, inclusion of the 
curvature will modify scaling functions Fj and G/, which will now depend on two variables: Fj{2A/ (u^uj), A/ Ac) 
and Gi{2A/ (u^Lu), A/Ac). We have not attempted to determine the most general form of these functions. However, 
we can make certain conclusions about their behavior near the mass shell. Indeed, as power-law singularities are 
cut at Ac, a particular contribution to the self-energy for A <^ Ac is obtained by taking an explicit result for this 
contribution for A Ac and replacing A by Ac. For example, the particle- hole contribution, given by Eq. (|3.29|) for 
A Ac, takes the following form for A <C Ac: 



ImS^^H — 



47r Ef 



(3.46) 
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Function Gi{x) in this form is still given by Eq. H3.30|l . We recall that Gi{x) behaves as Inla;!"""^ and a;ln|a;| for 
\x\ 3> 1 and |a;| <C 1, correspondingly. Using the small- a; asymptotic form of Gj{x), we find that the second term in 
Eq. H3.46|l is much smaller than the first one for |w| <C oJc- Therefore, ImSp^ in this limit is given by 



ImI]^H(^) ^ 111 -St ' for ^ « (3-47) 

47r hjp u^\uo\ 



The opposite limit of \uj\ ^ ojc (large x) exists only for a finite curvature. In this limit, ImEpg reduces to 

ImS^H(w) = 7^ -F^ In for ujli>uj,. (3.48) 

27r tiF |w| 

A similar procedure is applied to the contribution from ImE^. Away from the mass-shell, ImE^ is given by Eq. H3.25|l . 
A finite curvature cuts the infrared logarithmic divergence in the same way as in the second-order diagram. As a 
result, we obtain on the mass shell 

To logarithmic accuracy, a general form of ImS^ can be written as 

ImSf,(a;, A) = - In , ^ , . (3.50) 

''''^ ' %TT Ep max{|A|,Ac} 

Finally, the scaling function for the zero-sound contribution [Fi{x) from Eq. H3.36|l ] is small as a power-law of either 
X (for small a;) or x~^ (for large x). Therefore, for both \uj\ lOc and \uj\ 3> ujc regimes, the zero-sound contribution 
ImE^s can be neglected compared to ImEpjj + ImE^. 

Combining the formulas for the mass-shell forms of ImS^jj [Eqs. H3.47|) and (|3.48|l ] and ImS^ [ Eq. (|3.49|) ]. we 
arrive at 



up' r llnw^l, for |w|<Wc; 
^ ~ Ef^ \ \-i^EFl\'jj\, for |w| > 



I-S^ - X ,J - - ' - ' , — (3.51) 



Notice that for |w| <C Wc, there is no w- dependence in the logarithm, i.e., the frequency dependence of ImSp is 
perfectly regular in this range of w. Notice also that IniEp (w) remains positive for all frequencies, as it should in 
order for the quasi-particles to be stable. 



5. final result for imaginary part of the self-energy on the mass shell 



a. contact potential The net self-energy is a sum of forward scattering and backscattering contributions. On 
the mass shell (A = 0), forward- and backscattering contributions to ImS^ are equal to each other for ^ ujc [see 
Eqs. (|3.5|) and H3.51|l ]. whereas in the opposite limit of \uj\ <C ojc, the forward-scattering part [ Eq. H3.51|l ] is smaller 
by a large logarithmic factor than the backscattering one. The leading-order result for the on-shell ImE^ can then 
be written as 

In.E- (.) . ,„,E« M + I,nE« (.) ^ £|i 1„ {-^) , (3.52, 

where 

As we see, the non-perturbative effect in ImE^ on the mass shell is rather benign: all we have is a smooth crossover 
function interpolating between two different values of the numerical prefactor in a familiar In |w|-dependence 28]. 
Away from the mass shell (at A ^ 0), the non-perturbative effect in ImS^ is much more pronounced, and a non- 
monotonic behavior of the zero-sound term (|3.35l rOH|l gives rise to a non-monotonic variation of ImE^ near A = 0. 
We will return to this issue in Sec. IIII Fl where we discuss the spectral function. 
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b. finite range potential For a finite-range potential, a factor of in the backscattering contribution is [l^ 
replaced by u§ + u^i^^ — UQU2kF, where 

uo = mU (0) /2tt, U2kp = mf/(2fci.)/27r. (3.54) 

In the forward-scattering contribution, u is just replaced by uq. As a result, the net imaginary part of the self-energy 
on the mass shell becomes 

ImS« (c.) = ImSf (c.) + ImSf (-) = ^ In (^^) , (3-55) 

where -Umax = max{uo, U2feF}, and 

i/o\%°^w''"^r"""°^\^T'''t/^^ (3-56) 

[l/2-h(2Mo) 'M2feF(u2feF - Uo), for x < 1. ^ ^ 

E. Real part of the self-energy 

Next, we consider what happens to Rel]^(tj) near the mass shell. For definiteness, we consider a; > but A = uj — ek 
can be of any sign. The real part of the self-energy can be obtained either by a Kramers-Kronig transformation of ImE^ 
or directly, by evaluating the self-energy in Matsubara frequencies and analytically continuing it to real frequencies. 

1. backscattering 

First, we present the result for the total backscattering contribution to the self-energy. (By "total", we mean the 
sum of g2~ and 2kp— contributions). A Kramers-Kronig transformation of Eq. H3.5() yields, on the mass shell, 



ReEg (u;) = -V / ^^j^ '- = (3.57) 

n J h — u! 8 hp 

A non-analytic, cjlwj-behavior of ReEg (ui) is obviously related to a non-analytic, oj'^ In |w|-behavior of ImE3 (u;, k). 

2. forward scattering 

The real part of the self-energy consists of three contributions: from the remainder term (S^), from the particle-hole 
continuum (Sp^), and the from the collective mode (Sf^g). 

a. remainder Performing Kramers-Kronig transformation of ImS^ (given by H3.26|l 'l. we find 

Rel]fjw,fc) = -——V / dz^— In- ^ 



-V hm rd.i:i±^ln.^- 



Stt^Ef Ef^ooJ_^^ x~A \x\ 

The first term in Eq. 1)3.5811 is responsible for the renormalization of the effective mass and Z— factor, and we neglect 
it. The second term has the right frequency dependence, but its value on the mass shell depends on how we 

take the limit A ^ 0. This ambiguity is due to the logarithmic singularity in ImE^ [cf. Eq. H3.26|l ]. To eliminate this 
ambiguity, one has to re-evaluate the integral using the full form of ImS^, obtained by keeping the curvature finite. 
This form is given by Eq. I|3.49|) and is independent of A for A ^ 0. Performing a Kramers-Kronig transformation of 
Eq. I|3.49|l . we obtain 

ReSf.= ^^^ + 0(A2logA). (3.59) 
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Alternatively, one could just notice that 



In \iu\ — Im I — In [—{u 



Z0+)21 



Substituting Eq. H:i60(l into Eq. (|3.49() . one obtains the fuh (complex) S^: 



Stt Ef 



El 



(t^ + i0+)2 



(3.60) 



(3.61) 



Taking the real part of Eq. 13.61|l . we indeed reproduce Eq. (|3.59() . 

b. particle-hole contribution The same reasoning can be applied to the particle-hole contribution, Spjj. The 
imaginary part of Spjj near the mass shell is given by Eq. (|3.33|) . Using relation (|3.60f) again, we restore the full Sp^ 
as 



= i In 

8n Ef 



El 



. 1^1 A , 
+ 1— — In 

AttEf 



A2 



where we have also kept a first sub-leading term in A. The real part of Eq. H3.62|l is given by 

lu\lu\ |a;|A 



8 Ep AEp 



(3.62) 



(3.63) 



Adding up Eqs. H3.59|) and H3.63|l . we see that the leading oj |cj|-terms cancel each other, whereas the rest vanishes 
linearly on the mass shell: 



ReE^ 



PH 



ReS, 



MA 

4£;f ■ 



(3.64) 



Absence of a non-analytic, cjjcijl-term in ReEpjj + ReE^ is consistent with our earlier observation that on the mass 
shell ImEpjj -l-ImE^ is an analytic function of frequency [it scales as uj^v?' Inu, see (|3.51ll ]. Notice also that Eq. (|3.64|) 
is independent of u . On its own, such a term in the self-energy will give rise to the linear-in-w and u-independent 
correction to the density of states. We will see, however, that this term will be cancelled out by the contribution from 
the zero-sound collective mode, so that the full density of states remains analytic in uj. 

c. zero-sound contribution Next, we consider the contribution from the zero-sound collective mode. The real 
part of S^g can be obtained either by a Kramers-Kronig transformation of Eq. (|3.35|) . or directly from Eq. Ij3.23b|l . 
by expanding 1 — UIl{q) near the pole and performing the frequency and angular integrations. Either way, we obtain 
for A < w 



ReSfs = ^Re 



QdQ 



- 2A/(wFg) 



Evaluating the integral, we obtain 



ReE« 



zs — "S" ~Er- Er , 

8 Ef Xw'uj 



2A 



where 



Fr(x) = Re 



1 + VT^' 



(1 H — x)y/l — X -\ — x" In ■ 

2 2 \/—x 



(3.65) 



(3.66) 



(3.67) 



Subindex R implies that this is the scaling function for ReE^g- A plot of Fr{x) is presented in Fig. [S] Contrary 
to Fi(x) [ Eq. H3.36|l ]. scaling function Fji(x) is not one-sided, i.e., it is nonzero for both positive and negative x. 
However, it is clear from the plot that this function is asymmetric with respect to x. In fact, Fr{x) = for x > 1. 
For large and negative x, 



Fr{x) « 4Re[l/(5V^)], 



(3.68) 



which means that 



ReEfs oc u^Lo'^ {uj/ 



a/2 



(3.69) 
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FIG. 5: Scaling function Fr (x) [ Eq. JHUZJ]. 

for A < 0. This is consistent with the large x behavior of Fi{x) ~ Im[27r/ (5\/— x)]. Obviously, for large x, full S^g 
can be written as 



^zs — 



U^^/2 Uj\uj\ ( CJ\l/2 



20 Ef 



A 



(3.70) 



We see that away from the mass shell, i.e., at |A| ^ uj, the collective-mode component of ReE-'^(fc, w) scales as 
M'^w^, which is smaller by a factor of u compared to the u^w^ -contribution to the self-energy from backscattering 
( Eq. H3.57|l ). This smallness is another consequence of the fact that quanta of zero sound are not free bosons: the 
residue of the corresponding propagator scales as u^w^ and is thus small. For x — s- 0, 



Fr (x) 1 + a; + -a;^ In x ^ + 

On the mass shell, i.e., for x = 0, the function Fr approaches a finite value Fr = 1, so that 



ReE^ 



ZSlA=0 



v? ujyL)\ 



(3.71) 



(3.72) 



We see that the real part of the self-energy due to the interaction with the collective mode is strongly enhanced near 
the mass shell, such that at A = one power of the small parameter u is eliminated, and ReS^g becomes of the same 
order as the self- energy due to backscattering. This is one of the central results of this paper. Keeping the linear- in- a; 
term in Eq. (|3.71() results in a linear-in-A correction to the self-energy 



ReS^ 



uj\uj\ |w|A 



zs 



(3.73) 



8 AEp 

Adding up this result with Eq. H3.64|l . we find for the total contribution to the on-shell ReS^ from forward scattering 



ReEf 



ReEpjj 



ReEf 



ReE^s = - 



V? uj\uj\ 



Ef 



(3.74) 



3. final result for the real part of the self-energy 



a. contact-potential Combining the backscattering and forward scattering contributions to the self-energy for 
the contact potential [Eqs. H3.57|l and H3.74|l . respectively], we find that non-analytic (w and wA) terms cancel out. 
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and the net self-energy vanishes on the mass shelh 

RcS^ = ReEf + ReE^ = O {v? In A) . (3.75) 

This is another central result of the paper. It means that a non-perturbative contribution of the zero-sound mode 
totally changes the result of the second order perturbation theory, where to order v? we had ReE''^ oc iP'lo . 
b. finite-range potential For a finite-range potential, the backscattering part of the self-energy changes to 



RcEf = -|M {ul + ul,^-uoU2k,), (3.76) 

where uq and U2kF are given by Eq. (|3.54|l . The forward-scattering contribution comes only with Uq and is obtained 
from Eq. 13.74|l by replacing u ^ uq. A cancellation between backward and forward-scattering parts of ReE^ is no 
longer in place, and the net ReE^ is given by 

ReE« = |M U2fc, (7/o-"2fcJ. (3.77) 

For Uq 7^ U2kF , it is a non-analytic function of uj on the mass shell. 

Finally, in the ZS contributions to the self-energy. Eg. H3. 35(1 . u is replaced by uq as this contribution comes only 
from forward scattering. 

F. Spectral function 

1. short-range potential 

A non-monotonic variation in the zero-sound part of the self-energy is manifested in a specific feature in the spectral 
function 



A{iu,k) = --IniG^(w,fc) 

TT 



ImE^ {uj, k) 



TT [A -f ReE« (w, k)f + [ImE« (w, k)f ' 



(3.78) 



Having in mind a potential comparison with the experiment, we present a detailed discussion of A{uj, k) in this Section. 

A variation of ReE^ (w, k) has only a little effect on the shape of the spectral function, and we verified that it can 
be safely ignored. The effect of ImE^ {lu, k) is much stronger. As we have shown in Sec. IIII Dl ImE^ (cu, k) is a sum 
of four contributions 

ImE^ {lu, k) = ImEf + ImEf^ + ImEfn + ImEfg. (3.79) 

The particle-hole and zero-sound contributions contain scaling functions G/ and Fj which evolve as a function of 
A = w — efc on a scale A ~ A* = v?'oj/2 <C oj. The backscattering part, on the other hand, evolves only on much 
larger scale: A ~ w (cf. Appendix A). Therefore, one can safely put A = in ImEg , i.e., use its mass-shell value given 
by Eq. 13. 5|) . Finally, ImE^ crosses over between the forms given by Eqs. (I3.26f) and (|3.49ll at A ~ Ac = /Ep. 
Two situations are then possible, depending on the relation between Ac and A*. If Ac ^ A* or, equivalently, 
w ^ (jjc = v?Ep, the variation of ImE^g and ImEpjj occurs in the range where ImE^ can be approximated by 
its small- A form [ Eq. H3.49|l ]. which is independent of A. The sum of ImE^ and ImEg then vanishes, so that 
ImE^(cij, k) = ImEpjj + ImE^g. This sum can be further decomposed as 

ImE^(w,fc) Fo + ri(a;). (3.80) 

where x = 2A/u^w = A/ A*, Fq is a A— independent part of the self-energy [the first term in the particle-hole 
contribution, Eq. H3.29|l ]: 

To-f^ln^, (3.81) 
47r hp M l^^l 
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FIG. 6: Main panel: a log-plot of the spectral function A{io,k) [ Eq. (I3.83ll 1 in units of 1/-r'^v?Ef as a function of a; = 
2[u> — (.k)/v?Lu for — 2 and 7 — 0.05. A kink at a: = 1 is due to the interaction of fermions with the zero-sound mode. Inset: 
part of the spectral function A\{u,k) [ Eq. I3.85cll 1 in the same units function of x for 7 = 0.25. A maximum in yli at x = 1 
gives rise to a kink in total A (main panel). 



and ri(a;) is a sum of the scaling terms in ImSpjj and ImSf^g [Eqs. (|3.29|) and H3.35|l . respectively] 

T,{x) = ^^[Gi{x)+Fi{x)]. (3.82) 

Notice that F/(0) = and 0/(0) = 0. Substituting Eq. Ipn^ - lpTl^ into Eq. itXT^ and neglecting ReS^, we obtain 
a scaling form of the spectral function 

A( i.\ 1 L^ + Gi{x)+Fjix) 

A[uj, k) = „ „ ^ 5", (3.83) 

n^u^EFx^+-f^L^ + Gi{x)+Fiix)f 

where 

and = In (^Ep /u^\u;\') is a large factor. We consider a setup when oj is fixed and the spectral function is measured as 
a function of the momentum. This is equivalent to varying x at fixed lu in Eq. H3.83|l . In photoemission measurements, 
this setup produces what is known as a "momentum distribution curve" (MDC). A plot of A [ Eq. (|3.83|) ] as a function 
of X is shown in the main panel of Fig. For solely illustrative purposes, we have chosen — 2 and 7 = 0.05. 
We see that the spectral function contains not only a narrow quasi-particle peak at a; = (i.e, uj = Ck) but also a 
well-pronounced kink at a; 1. To understand the reasons for the kink in A{lo, k), we notice that at typical A ~ A*, 
Fi is of order uP' / Ep, which is smaller than Fq by a large L^^. Therefore, we can expand the spectral function in Fi 
and represent ^(w, k) as a sum of two contributions 



A{oj,k) = Ao(w,fc) -f Ai(w,fc) (3.85a) 

1 U 

tt'^vPEp x^ -I- ^\ 



M^.k) = , " , (3.85b) 



Ai(^,fc) = -^1— [G,(x)+fj(x)] /:' (3.85c) 

where 71 = ^L^. The first term Eq. (|3.85b|) describes a regular quasi-particle peak at A = of width 71. The scaling 
behavior of the self-energy shows up in the second term, Eq. (|3.85cf) . Since 71 <C 1, damping affects the behavior 
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of Ai {uj,k) only at very small x: cc ~ 71 ^ 1. For these a;, Ai exhibits a rapid non-monotonic variation, but it is 
overshadowed by the rapid variation of Aq. For |x| ^71, ^0 is smooth, whereas Ai is determined by the sum of the 
scaling functions, which is non-monotonic in x 



Ai {uj,k) 



1 



Gi (x) + Fi (x) 



u^Ef 



(3.86) 



Comparing now the behavior of two scaling functions, G/ and F/, we see from Fig. 0] that Gi{x) varies smoothly at 
a; ~ 1, whereas Fi{x) has a sharp peak at a; = 1. We remind (cf. discussion in Sec. IIII D)| that this sharp peak is 
associated with the fact that, for |A| > A*, a fermion with energy w > above the Fermi level can emit ZS bosons 
with any frequency in the interval < f2 < w (or absorb bosons in the interval < 57 < — w for w < 0), whereas 
for |A| < A*, a Cherenkov-type restriction makes it impossible to emit and absorb bosons with frequencies above 
|A|/(1 — vp/c) « 2|A|/u2. The sharp peak in Fi{x) gives rise to a peak in Ai {cu,k) at a; = 1, see inset in Fig. 
The peak in Ai gives rise to a kink in the full A{lu, k) aX x — 1. We emphasize that the kink originates from the 
zero-sound contribution to the self-energy, i.e., it reflects an essentially non-perturbative effect. 



'Ep. In this case, for x ~ 1, ImE^^ can be replaced 



The second situation occurs when A^, 3> A*, i.e., uj <^ujc 
by its large- A form [ Eq. H3.26|l ] . Re-expressing ImS^ in terms of the dimensionless variable x we obtain 



8^Ef 



In 



Ef 



In la; 



(3.87) 



Decomposing again ImE^ into x— independent and x— dependent parts, we obtain instead of Eqs. (|3.85al3.85c|l 

A{uj,k) = Aaiuj,k) + Ai{uj,k) (3.88a) 

(3.88b) 



Ao{uj,k) 
Ai {Lu,k) 



1 lnEl/u^\uj\^ 



2-k'^u'^Ef a;2 -f 7I 
1 



2'k^u^Ef 



Gi [x) + Fi {x) + - In |a;| 



9 2 



where 
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1^1 El 



AttEf u'^\uj\ 



3 • 



(3.88c) 



(3., 



The behavior of A{uj,k) in this case is a bit more involved than for Ac < A*; nevertheless, the a general structure 
is the same as before: the spectral function has both a quasi-particle peak at a; = and a kink at a; = 1, due to the 
contribution from the zero-sound mode. 



2. Coulomb potential 

A kink in the spectral function due to the interaction of fermions with the collective mode is not a special feature of 
the model with a short-range repulsion but a general phenomenon. To illustrate this point, we consider a 2D system 

1 /2 

with the Coulomb interaction, when the collective mode is a plasmon with dispersion noiQ) = [e^mvlQ] ' . An 
on-shell electron can emit plasmons only if its energy exceeds a certain critical value: \u}\ — \ek\ > t^pi = V^rgEp 
[33 , where is the usual ideal-gas parameter for a charged system, which is assumed to be small. Although formally 
there is an interval of energies in between Wpi and Ef, in practice it cannot be very large. In what follows, we will 
consider only the low-energy limit: \lu\ ,\£k\ ^ i-^pi- In this case, emission of plasmons by electrons is possible only 
away from the mass shell (w 7^ ek), and the effect we are interested in is a kink in the spectral function rather than 
the lifetime of an electron. 

Near the plasmon pole, the Coulomb potential reduces to 

^(0,0). /";""^ (3.90) 

{n + id) - (Q) 

and, correspondingly, the imaginary part of the self-energy is given by {uj > 0) 

ImS^ (w, Ck) = 2Tre^mvj, J SQS {uj - ek - n - vfQ cos 0) 6 {n^ - nl (Q)) . (3.91) 
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The second S— function forces the boson momentum Q to be small: Q = /me^vj^ ~ (|r2|/vi?) /wpi) <C \^\ jvp- 
Therefore, one can neglect the Q-dependcnt term in the argument of the first 8— function. Performing an elementary 
integration and considering the case of w < in the same way, we obtain for ImS^ (w, efc) : 

\mY.^(tjj€ ) = [ ^'^k - i^f /Ep, min{0,w} < efc < max{0,a;}; 
' \ 0, otherwise. 

For fixed w, ImS^ has a kink at efe = 0, i.e., at the Fermi surface, where ImE^ vanishes discontinuously (see comment 
1^). (Keeping the Q-dependence in the first i5-function in Eq. (|3.91|l . one sees that in fact the kink and zero of 
ImE^ are separated by a small energy scale w^/wpj.) Notice that the electron charge dropped out of the result. 
It can be shown that the spectral function of a bi-layer system with two plasmon modes-with ^JQ and acoustic 
dispersions-behaves in a similar way but we defer a detailed discussion of this case to later occassion. 

A kink in the spectral function could, in principle, be detected in photoemission experiments on layered materials 
[29| or in a momentum-conserving tunneling between two parallel layers of 2D gases ^] . 



3. absence of a non-analytic correction to the tunneling density of states 



Both the particle-hole and zero-sound contributions to the real part of the self-energy contain a specific term, which 
is proportional to the product |a;| A and is independent of the interaction [cf. Eqs. H3.63|l and (|3.73() ]. Each of these 
terms on its own would give rise to a linear-in-|a;| and u-indepcndcnt correction to the density of states. Indeed, a 
term in the self-energy of the form 

s\uj\/:\/Ef, (3.92) 
where s is a numerical coefficient, gives rise to a non-analytic frequency dependence of the renormalization factor: 



Z{uj)=\l + s^j (3.93) 

[We defined Z{u)) in such a way that G^{k,uj) = Z{uj)/{uj - Ck + iO+).] The linear -in-jcijl-term in Z{ll!) results in a 
linear frequency dependence of the tunneling density of states 

Niu) ^-^ f -f^IniG^ [u., fc) = - Ziu) « {mM (l - sM.) . (3.94) 
Ti" J (27r) Ti" \ Ef J 

However, we see that the A|(jj|-terms in Eqs. (|3.t)4|) and (|3.73|) cancel out in full ReS-'^ [ Eq. (|3.74|) ]. so that s = 0. 
Therefore, to order m^, N{lu) is analytic in uj. This result is valid for any finite-range potential as the cancellation of 
I a; I A— terms occurs between the forward-scattering contributions to the self-energy, all of which contain the same 
coupling uq. Our result that there is no linear-in correction to the DOS is in agreement with Ref. 39|, where the 
tunneling density of states was obtained for the case of a multi-layer system with the Coulomb interaction. Inter-layer 
screening gives rise to an acoustic branch of the plasmon spectrum which is an analog of the ZS mode of our model. 
Notice that an |a;|-correction to the density of states does exist for a single layer with the Coulomb potential |39l.l40l|. 



IV. SPECIFIC HEAT 



As we pointed out in the Introduction, the a;|t£j-non-analyticity in the second-order self-energy gives rise to a non- 
analytic, u^T-correction to the ratio C (T) /T [l^. However, it was shown in Sec. IIIIEl that. upon re-summation, 
higher-order forward-scattering contributions to ReS^ also becomes of the order near the mass shell and modify 
the second order result. For a contact interaction, the non-perturbative contribution even cancels the second-order 
M^a;|aj|-term in ReE^. The question addressed in this Section is whether the forward-scattering component of the 
self-energy modifies the non- analytic term in the specific heat. We show in several ways that this does not happen, 
i.e., the enhancement of the forward scattering self-energy near the mass shell does not affect the specific heat. 

We also go beyond the weak-coupling limit in this Section, and consider the non-analytic behavior of the specific 
heat in a generic Fermi liquid. 
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A. Specific lieat via self-energy 



A relation between the entropy (and thus the specific heat) and an exact fermion Green's function can be found 
in Ref. P]. However, this relation is justified only for the Fermi- liquid, linear- in-T part of C(T). In order to find a 
sub- leading, non-analytic contribution to C(T), one needs to re-examine the assumptions, made in Ref. 0, and to 
establish a correct relation between C(T) and the self-energy beyond the leading order in T. 

The relation between the thermodynamic potential and the Green's function reads 

S = 2r^ /'^lnG(a;„,fc,T = 0), (4.1) 

where G{ujm,k,T = 0) is the Green's function evaluated at discrete Matsubara frequencies but with no additional 
T-dependence. Converting the Matsubara sum into the contour integral and using a familiar thermodynamic relation 



C{T) = -T—, 



one obtains 



C(T)/T = - 



2T d 
~ df 



(27r)2 



dujuj^^ arg G^(uj, k) 

OU! 



(4.2) 



(4.3) 



(Notice that there is no need to distinguish between Cp and Cy here, as we are interested in the T^— term in 
C (r), whereas the difference between Cp and Cv is of the order |^.) Both the retarded and advanced Green's 
functions in Eq. (|4.3|l are evaluated at T = 0, thus the derivative in the r.h.s. of Eq. H4.3|l affects only the Fermi 
function [the temperature derivative of ng was converted into the frequency dependence by using a familiar identity: 
dno/dT = -{uj/T)dno/duj]. 

Eq. H4.3|l correctly describes the regular, linear-in-T part of the specific heat. Indeed, substituting the analytic, FL 
form of the self-energy, Ea. p.l|l . into Eq. 1)4. 3|l . one finds that C(T) is given by the Fermi-gas result [ Eq. but 
with a renormalized mass, which is composed from coefficients a and b in Eq. 

It was conjectured in Ref. 1] that Eq. (|4.3I) describes not only the leading but also the sub- leading terms in C{T). 
However, this conjecture is questionable, as the accuracy of the low-temperature expansion used in the derivation of 
Eq. (|4.3|) is not specified. In other words, it is not obvious that if one retains contain higher powers of oj in S, one 
should not at the same time take into account an explicit temperature dependence of S. Of particular concern are 
the situations when the self-energy depends on T via a scaling function of variable uj/T (this happens in our case; 
see below). As typical oj are of order T, the argument of the scaling function is of order unity, thus neglecting the T- 
dependence is not justified. Moreover, it was shown in Ref.*^ that the zero-temperature and temperature-dependent 
parts of the self-energy give comparable contributions to the non-analytic, InT part of C{T) in 3D. 

We will still be considering the case of a weak interaction, when |E^| <C |a;|. In this case, Eq. H4.3|l can be simplified 
further by expanding the logs of Green's function in the self-energy, which results in 



C{T)^Cyg{T)+6C{T), 
where C-pg{T) is the specific heat for free fermions in 2D [ Eq. Hl.ll|l ] and SC{T) is given by 



2 d 




r (fk r 


n df 


T J 


' (2^)2 1 



(4.4) 



(4.5) 



There are two contributions to 5C{T)-one from ReE^ and another from ImS-'^-which we label as Ci(T) and C2(T), 
correspondingly: 



SC{T) = Ci{T) + C2{T) 
d 



C,{T)/T = -2 ^ 
2 8 

C2{T)/T = - ^ 



(Pk 

(2^ 



diouj^V ^^Iml]^(w, fc, T = 0) 



(4.6a) 
(4.6b) 

(4.6c) 



In the expression for Ci(T), we have used the fact that lmGQ-{uj — k) = —tt5{u! — e^). The delta- function in l|4.6b|) 
implies that Ci(T) is determined by ReS-'^ only on the mass shell, where cu = ek- The second term C2(T) contains 
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the integral of lmT,^{u!, k). If lmYi^{uj,k) depends on lu but not k, which is the case, e.g., for the elect ron-phonon 
interaction 41], the momentum integral in Eq. (|4.6cl) vanishes once one approximates the density of states by a 
constant, so that C^iT) drops out. Indeed, 



C2 (T) cx / dL0Uj^lm5:^{uj)V 



dek- 



w - efc 



= 0. 



(4.7) 



However, for a general case, when S depends on both uj and fc, there are no a priori reasons for C2(T) to vanish, and 
thus the imaginary part of the self-energy contributes to the specific heat as well. In what follows, we will omit the 
analytic terms in which just renormalizes the coefficient of the linear T-dependence in Eq. Hl.ll|l . and consider 
only the non-analytic contributions to Ci and C2. 

Next, we compare the result for the correction to the specific heat given by Eqs. (|4.6al4.6c|) to the one obtained in 
a different way, namely, employing the Luttinger-Ward formula for the thermodynamic potential S: 



d^k 



ln(GoG"^) 



(4.8) 



Here, .^0 is the thermodynamic potential of the free Fermi gas per unit area, Gq — {iuj„i ~ ^k) i G ~ 
{iojm — Cfc + ^, S is the exact (to all orders in the interaction) self-energy, and is the skeleton self-energy 
of order i/. Both the skeleton and full self-energy, related via 



(4.9) 



are evaluated at finite T. The diagrams for J^^, are obtained from those in Fig. El by replacing the bare Green's 
function and interaction lines by the exact ones. Expanding both G and in Eq. H4.8|l back in E, one generates a 
perturbative expansion for S. To second order, diagrams generated by the first two terms in Eq. 14.8|l correspond to 
self-energy insertions into a free thermodynamic potential (circle). Such diagrams just renormalize the prefactor of 
the leading, linear-in-T part of C{T). The non-analytic contributions come from the third (skeleton) term. To second 
order in the interaction, this contribution is 



SE 



fk 
4^ 



T,{uJrn,k,T)Go{k,UJ„ 



(4.10) 



where T,(LOrm k,T) is (a non-analytic part of) the second-order self-energy. Converting the Matsubara sum to an 
integral over real frequencies and using relation Eq. H4.2|) between S and C{T), one obtains 



d^k 
(2^ 



duj Im{G^(cj,fc)E^(w,fc,r)} 



no{uj) ~ - 



(4.11) 



We emphasize that in this approach E^(fc, cjm, T) is evaluated at finite temperature. 

Generally speaking, Eq. H4.11|l and Eqs. H4.6al4.6c|l give different results for SC. Indeed, let us assume for a moment 
that E^ depends only on frequency but not on k and T. 



Then 5C{T) from Eq. l|46H|l and from Eq. differ by 

a factor of four. The derivation based on the Luttinger-Ward formula is free from assumptions on what constitutes 
the main source of the T-dependence in 5C(T). In fact, Eq. (|4.11() is valid for any temperature albeit for weak 
interactions. The safe way to proceed therefore is to use Eq. (j4.11|l but not Eq. 1)4. 5|l . It appears, though that for our 
case, there exists a deeper relation between the two formulas. Namely, the two expressions yield identical results for 
SC{T), provided that one replaces E(cj, k,T — 0) by the temperature-dependent self-energy in Eq. I|4.5|) . i.e. 



m^6CiT)/T=- ^ 



(4J2) 



This is how SC{T) was calculated in Ref. j^. The overall factor of four difference between Eq. (I4.11|) and Eq. 14. 5|) 
is eliminated by two reasons. First, the imaginary part of the self-energy to order does depend on ek albeit only 
logarithmically: ImE^(a;, k) ex. uP' In \lo + tk\- Then Eq. 14.6c|) gives the same contribution as Eq. (|4.6b|) . The details 
of this calculation are presented in Appendix |0 An additional factor of two appears because the derivative over T in 
Eq. (|4.12|) now acts not only on tiq but also on E^. In Appendix lO we show that these two terms contribute equally 
to SC{T); hence, an additional factor of two. 
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To summarize, Eq. (|4.11() gives a correct result for a non-analytic term in C(T) to second order in the interaction 
without any assumptions or constraints. At the same time, Eq. H4.5|l gives the correct result provided that the self- 
energy in Eq. H4.5|l is evaluated at finite T rather than at T = 0, as specified by Eq. H4.12(l . This is a consequence 
of the uj/T scaling in the non-analytic part of S(ti;,T). We did not study, however, whether or not this statement is 
specific to our weak-coupling case or has a wider range of applicability. Having this precaution in mind, we will be 
using Eq. (|4.12(l in the following analysis. 

It is convenient now to separate the self-energy into the zero-sound part and the rest, which includes the 
backscattering- and PH-contributions from spin- and charge channels, as well as the remainder term, Eex- Such 
a separation is convenient because the imaginary part of ImS^ has a substantial fc-dependence and thus, according 
to the discussion in the previous Section, gives a contribution to the specific heat. On the other hand, the imaginary 
part of the rest of the self-energy depends on k only logarithmically, and will be shown not to contribute to C (T) . 



1. non- zero- sound contribution to C{T) 



We begin with the part of the self-energy that contains all contributions but the zero-sound one: 



■^PH 



Consider first the contribution to the specific heat from the real part of the self-energy, Ci(T), Eq. Ij4.6b|l . A sum 
of the two contributions, ReSpH + ReEcx, vanishes on the mass shell and therefore does not contribute to Ci(T). 
A non-analytic part of ReEg on the mass shell and at T = is given by Eq. H3.57|l . At finite temperatures, ReSg 
has been calculated in ^^l; the result of this calculation is that ReE3(fc,cj,T) differs from RcSg(fc, Li;,T = 0) by a 
multiplicative factor which is a scaling function of lo/T: 



ReEf(fc,w,T) = ReEf (fc,w,T = 0).g(tj/r); 



g{x) = 1 + — 

where Li2(x) is a polylogarithmic function. 

Substituting Eq. (|4.13|l into Eq. Ij4.6b|) we obtain [3 



h Ll2 

12 



(4.13) 



C^[T)IT = -^-^u^CfgIEf, 



(4.14) 



where Cpc is given by Eq. (|1.11|) . As is expected, the non-analytic u) |ci;|-dependence of ReEg gives rise to a non- 
analytic contribution to the specific heat Ci (T) /T cx v?T. 

To logarithmic accuracy, ImE(a;, k) oc In |ti;| does not depend on fc, hence, according to Eq. 14.7II . C2(T) vanishes. 
To demonstrate unambiguously that C2 (T) vanishes, one has to go a bit deeper and analyze ImS^(fc,a;) beyond the 
logarithmic accuracy, focusing specifically on the momentum dependence under the logarithm. In Appendix [CI we 
show that the contributions to C2 (T) from ImE^I and ImEpH + ImEex cancel each other, i.e., there is indeed no 
contribution to the specific heat from ImE^. 



2. zero-sound mode contribution to C{T) 

We now use the results for ReE^g and ImE^g from Sec. IllII and calculate the contribution to the specific heat 
from the zero-sound mode. We show that the contributions from ReE^g and ImE^g cancel each other, ie., that 
the zero-sound mode does not contribute to C(T) despite the fact that ReE^g is enhanced near the mass shell. As 
we only need to prove the cancellation, we just use zero-temperature forms of E^g; as we explained in the previous 
Section, the difference between the results for the specific heat found using the zero- or finite-temperature forms of 
the self-energy, is just an overall numerical factor. 

Substituting ReEfg from Eq. (|3.72|l into Eq. I)4.6b|l . we obtain 

Ci{T)/T^-Nu' (4.15) 

TT V, 



where 

dxx^ 9 
cosh 



-r2- = oC(3). (4.16) 
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Substituting next Eq. H3.35|l for lmS^g{k,uj) into Eq. (|4.6c|l . we obtain 



where 



Z{T) = — u^V 
2m 



duJLJ 



cosh 



2T 



u - ek 



Im 



QdQ 



y/u'^ - 2{uj - ek)/{vFQ) 



Re-scahng the variables and replacing / cPk by (m27r) J dck, we obtain from Eq. H4.18|l 



where 



f°° dP 
J{cj)^vi — Im 



duJLJ 



cosh w/(2r) 



Q^'^dQ 
^u-^Q - P 



(4.17) 



(4.18) 



(4.19) 



(4.20) 



The integration region over momenta P and Q are defined by the following conditions: P > u^Q and uo > vfP/2. 
Evaluating the integrals in H4.20|l over this region, we find 



Jiu;) 



to 



2u \vf , 

Substituting this result into Eq. (|4.19() . and then into (|4.17|) . we obtain 



C2{T)/T = - 



d 



Snvj, dT 



1 

r2 



cosh^w/(2T) 



= u^N 



(4.21) 



(4.22) 



Comparing Eqs. (|4.15() and (|4.22l) . we find that these two contributions to the specific heat cancel each other, i.e., 
there is no non-analytic contribution to the specific heat from the zero-sound collective mode to second-order in the 
interaction, despite the non-perturbative enhancement of Sfk near the mass shell. An absence of the zero-sound 
contribution to the specific heat is another central result of the paper. 

The final result for a weak, contact interaction is then given just by the perturbative contribution, Eq. H4.14|) : 



Ci (T) /T = - 



Cfg J^F- 



(4.23) 



3. finite-range potential 

For a finite-range potential, the forward-scattering part of the self-energy involves only uq. Hence, all cancellations 
discussed in the preceding Sections are still in place. In particular, the sum of ReS^ and ReEpjj still vanishes on the 
mass shell and the contributions from real and imaginary parts of S^s C (T) still cancel each other. In addition, 
the momentum integral in Eq. (|4.t)c|l of each of the three terms ImEg , ImEpjj, and ImS^ still vanishes. As a result, 
the specific heat is again determined by the real part of the self-energy from backscattering. The self-energy due to 
backscattering for a generic U{Q) is given by Eq. (|3.76|) . hence SC (T) becomes 

6C (T) /T = - ^ {ul + ulk^ - uoU2k,) Cfg/Ef. (4.24) 

This expression is the final result for the non-analytic correction to the specific heat to order u^. For ug = U2kp = u, 
it reduces to the contact-potential result, Eq. (|4.14|) . 
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k + q 




p + q 



k + q 1 + q 




FIG. 7: Diagrams for the thermodynamic potential containing maximum number of particle-hole bubbles. For the Coulomb 
potential, diagrams 1(a), 2(b) and 3(c) represent ring diagram series to third order in the interaction 



B. Specific heat via the thermodynamic potential 

The calculation of the specific heat via the self-energy presented in the previous Section is quite involved, as it 
requires a detailed knowledge of S (w, k). To verify the main result of the previous Section-that there is no contribution 
to the specific heat from the collective mode-and to understand it from a different perspective, we employ an alternate 
approach. In particular, we obtain the specific heat by finding the thermodynamic potential, S, directly, and then 
using relation Eq. (|4.2() . In this approach, fermions are integrated out from the very beginning, and the intricate 
details of their self-energy are not important. Another advantage of working with the thermodynamic potential is 
that the entire calculation can be performed in Matsubara frequencies. 



Luttinger- Ward expansion 



To generate a perturbative expansion of S, we follow the Luttinger- Ward approach |42,|43j, in which 2 is expressed 
in terms of the exact Green's functions and the skeleton self-energies, as specified in Eq. (|4.8|l . To begin with, we 
consider the contact-interaction case. The diagrams for the self-energy describing the interaction of fermions with 
collective modes have been discussed in Sec. IIIII Expanding the result for the self-energy, Eq. (|3.22|) . back in powers 
of interaction U and substituting the resulting series for E and T,^ into Eq. 14.8|l . we generate the series for the 
thermodynamic potential. The non-trivial diagrams for S up to third order in U are shown in Fig. [7| Explicitly, 



2c/n„ 



5 



■ in(i + un,n) 



(4.25) 



where Il„i{q) is the polarization bubble in Matsubara frequencies. Again, we will need only the asymptotic forms 
of n„j((7) for Q near and near 2kp. Analytically continuing the corresponding retarded expressions [Eqs. H3.3|l and 
(|3.4() ] ] to Matsubara frequencies flm — lirmT, we obtain 



27r 



(4.26) 
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and 



1 - 



Q-2kF 

2kF 




Q-2kF 

2kF 



1/2 



2kFVF 



,Q«2/cF. 



(4.27) 



One observes immediately that the series in Eq. H4.25|l converges for small U because n„j Q) is regular for any 
Vim and Q [in contrast with Q) which is singular at the boundary of the particle-hole continuum 57^ = [vfQY]- 

Every order in U then gives a finite contribution to the thermodynamic potential which can be calculated separately 
from other orders. Therefore, there is no need for re-summation of the perturbation theory for S for a weak interaction. 
This tells us that the 0{U'^) term in the specific heat cannot have any non-perturbative contributions, i.e., to second 
order in [/, S is given just by 



" — "0 — "2 



-^(^7^„)^ 



(4.28) 



This explains the absence of collective- mode contribution to C{T), 
Sec. llV A2l 



which we have demonstrated explicitly in 



2. Evaluation of the thermodynamic potential in Matsubara frequencies 

Next, we show how the non-analytic correction to the specific heat emerges in the Matsubara formalism. The 
T^-term in C(T) comes from a non-analytic, T^-piece in S, and we will be searching for this term in Eq. 14.28|l . 
We first show that the Q~iiitegral of 11^ taken over momenta near Q = and Q = 2kF contains a non-analytic, 
i7^1n|f2m| part. Indeed, squaring the small-Q form of [ Eq. H4.26|l ] and substituting the result into Eq. Il4.28|) . 
we find 



2n 



{2nY 



dQQ 



m 



nl In- 



nl + {vFQ)^ (2^)3t;2_ ™ 



(4.29) 



Similarly, the square of the second term in the bubble near 2/cf [Eq. H4.27|l ]. yields another fl^ In |il„i|-term with the 
same prefactor, as from the region of small Q. To see this, we substitute the square of Eq. H4.27(l into Eq. 14.28|l . 
re-define the integration variable as x = (Q — 2fci?)/(2fci?), and retain only the square of the second term in 11^^: 



i 



QdQ „2 



2tt 



n^„(f^m,Q) 



(2^)^ 



dx 



\2kFVF J 



(4.30) 



The precise limits of the integration over x are not important. Expanding in frequency, we obtain 



Zn 



{2^y 



A2kF? 



dx [ X + \x\ 



ikFvl-lxl 



{2t:)^v 



3, ,2 ™ 
F 



(4.31) 
(4.32) 



Comparing Eq. H4.32|l and Eq. (|4.29|l . we see that the non- analytic contributions from Q = and 2kF are equal. 
Substituting the sum of the two contributions into Eq. (|4.28f) . we obtain 



n^. = r E ^m^^7^ = —S{M), 



(4.33) 



where 



M M-l 

S (M) EE V TO^ In — = -M (M + 1) (2M + 1) In M - V (m + i f In (m + 1) 
^-^ m 6 ^-^ 



m=0 



(4.34) 
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and M = [Ep/2ttT] ^ 1. The choice of Ep as an upper limit in Eq. (|4.33() is completely arbitrary; since are looking 
for a universal contribution to the thermodynamic potential, the choice of cutoff is not important. Next, we use the 
Euler-Maclaurin formula : 

E /M = / - 2 [/(^^) + /(O)] + E B^.L "^-J-J. Li, (4.35) 

m=l p=l ^ 

where Bk are the Bernoulli coefficients and Z*^"^ is the n-th derivative of /. Applying this formula to Eq. 14.34|l . we 
see that the derivatives /'^^'^"^^(M) for k> 2 form a series in 1/M for large Af, whereas /^^-' (A/) and /'^^''^^'(O) give 
A/— independent contributions for fc > 1. Combining the result of the Euler-Maclaurin expansion with the first term 
in Eq. (|4.34l) and taking the limit of Af ^ oo, we arrive at 

M 11 1 

V TO^ In — -M^ M -a + + . . . (4.36) 

^ m 9 12 360A// ^ ' 

The AT-independent term, a, is represented by the following series 

1^1. d^P- 



= 0.0304 . . . (4.37) 



9 ^ (2]5)! ^Prfa;2p- 
111 1 



9 12 360 7560 

Although we have not been able to prove this analytically, we observe that, to very high accuracy, 

" = S^- (4-38) 

The same constant is obtained when calculating the specific heat in real frequencies, when the Matsubara sums are 
converted into ingegrals (see Appendix (EJ . We will thus treat relation H4.38|l as an exact one. 

Terms of order Af^, 1/Af, etc. in S (M) generate regular-T", T^, T'^, etc. -corrections to S2, whereas a constant 
term (—a) gives a universal, non-analytic T'^— contribution to S2. This contribution is precisely what we need. 
Retaining only this term in Eq. H4.33() and substituting the result into Eq. (|4.2I) , we obtain the same correction to the 
specific heat as the one found by expressing the specific heat via the self-energy, Eq. (|4.23|) . 

For a finite-range interaction, a slight modification of the analysis presented in this Section leads to the result 
identical to that in Eq. (|4.24|) . 

For the sake of completeness, in Appendix ^ we evaluate the specific heat by computing the thermodynamic 
potential in real frequencies and show explicitly how the contribution from the zero-sound is cancelled out. 



C. specific heat in a generic Fermi liquid 

Now we are in a position to discuss a more general question-what happens to the term in the specific heat if 
the interaction is not weak. First, we discuss a model case of contact interaction of arbitrary strength and then move 
on to the case of a generic Fermi liquid. 



1. contact interaction 



To second order in contact interaction u, relevant diagrams for the thermodynamic potential contain the square of 
the polarization bubble [cf. Fig. [7| 2(a) and 2(b)]. Since Q — and Q = 2fci? contributions to the thermodynamic 
potential are identical for this case, we evaluate the Q = 0-contribution first, and then just double the result at the 
end. In Sec. IIV B 21 we have shown that the non-analytic, contribution to the thermodynamic potential comes 
from the square of the dynamic part of According to Eq. (|4.25|) . higher orders in u generate higher powers of n,„. 
Another effect of higher orders is that self-energy insertions result in replacing the bare Green's functions by exact 
one. As the main contribution to the T^-term in S comes from the states near the Fermi surface, one can approximate 
exact G by its expression near the pole 



G{uJm,k) = -, 



(4.39) 
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where Z is the renormaUzation factor, = v*p(k — kp), Vp = kp/m,*, and m* is the renormalized mass. Parameters 
Z and m* are some functions of the bare interaction u, whose forms, in general, are not known. This amounts to 
replacing the prefactor and the Fermi velocity in Eq. H4.26|l 



Tlm{nrn,Q)^Tl*^{n,Q) 



m*Z^ 
2n 



(4.40) 



A term of order n in series Eq. 14.25|l contains (H^)" and thus generates a binomial expansion in powers of V* 



(n;j" = i-nz^m*/27rrY.i~yciv*\ 



(4.41) 



1=0 



where is the binomial coefficient. It is easy to make sure that only the term with Z = 2 in Eq. (|4.41|) yields 
ri^lnjOml upon the momentum integration, whereas all other terms yield just fi^. As was shown in Sec. IIVB'2l the 

frequency sum 



S1,„=0 



(4.42) 



gives a universal contribution to the thermodynamic potential, responsible for the term in the specific heat 
(see H4.33|) '). At the same time, the frequency sum 



n,„=Q 



(4.43) 



contributes only analytic-T^, T^, etc. -terms to S, but no T^— term. The problem therefore reduces to collecting the 
combinatorial coefficients of V*^ terms at each order and re-summing the perturbation series. 
Expanding each term in Eq. 14.25|l to order V*^ and using 

k=l ^ 

we find that the T^- term in the specific heat is given by 

SC (T) /T : 

Here = 7rm*r/3, E*p ^ kFV*p/2, and 



9 "cff^FG/-^F- 



(4.44) 



(4.45) 



where 



u* = Z'^m*U/2TT. 



2 V 1 - u* 



2 V 1 + It* 



(4.46) 



(4.47) 



This result is valid for < m* < 1. The divergence of u^^ at u* = 1, resulting from the spin-channel, signals an 
instability towards a magnetically-ordered state. A plot of u^g (m*) is presented in Fig. |S1 



2. generic interaction 



The result for a contact interaction, Eq. 14.45(1 . is of a rather limited use, as, in general, the Fourier transform of the 
interaction does depend on the momentum transfer. To obtain a form of (5C(T), valid for a generic Fermi liquid, we 
first go back to the second-order diagrams for the thermodynamic potential, and identify the structure of the vertices 
contributing to the non-analytic part of C{T). 
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FIG. 8: Effective coupling for a T^-term in the specific heat u^g [ Eq. H4.46^ ] as a function of the renormahzed interaction u* , 
Eq. llOTji . 



We consider first a small-Q contribution to the diagram 2(b) in Fig. [7| It is proportional to 
U^{Q) TJ2nJ d^Q^^{^m,Q), where the integration is restricted to small Q. Each of the two polarization bub- 
bles is obtained by the integration over internal fermion momenta, k in the upper bubble and p in the lower one. At 
first glance, k and p are completely uncorrelated, as the integrations over k and p are independent of each other. 
This would imply that the total momenta for the two vertices in diagram 2(b) are arbitrary. In general, this is indeed 
true. However, a non-analytic, T'^-term in the thermodynamic potential arises only from a product of non-analytic, 
\^m\/Q parts of the two polarization bubbles, integrated over Q. In Appendix IfI it is shown that the |f2m|/Q-term 
in n,„(r2,„, Q) comes from an integral / (PkG {ujm + ^m, k -I- Q) G{u!m, k) from only over those regions of k where 
k is nearly orthogonal to Q (and the same for p in the other bubble). Since both k and p are almost orthogonal 
to Q, they are either nearly parallel or nearly antiparallel to each other. In Appendix^ it is also shown that the 
contribution from the nearly parallel k and p vanishes, i.e., the non-analytic , "Q — 0" contribution to S involves only 
a vertex with a small momentum transfer Q and small total momentum k -f p ( "backscattering" vertex) . The same 
consideration holds for the J7^(2/cF)-term from diagram 2(b) in Fig. 13 and for the C/(0)C/(2/cF)-term from diagram 
2(a) in Fig. In both cases, the T'^-term in S comes only from those momentum range, where k and p are nearly 
antiparallel. We thus see that the T^-term in C{T) involves only vertices of the type (k, — k; k, — k) and (k, — k; — k, k). 

Consider now what happens when we add higher-order terms in the interaction. They lead to two types of correc- 
tions: self-energy insertions into the fermion lines in the two bubbles and corrections to the vertices. The self-energy 
corrections are of the Fermi-liquid type: they account for the appearance of the quasi-particle Z-factors, and for the 
replacement of the bare fermion mass by the effective one. Vertex corrections lead to a variety of diagrams. A typical 
n-th order diagram contains n bubbles. [We remind that these bubbles are either explicit, as in diagrams 2(b) and 
3(c), or are obtained after integrating over the fermion variables, as in the rest of the diagrams. Some diagrams, e.g., 
3(b), contains bubbles both explicitly and implicitly.] To obtain a contribution to C{T), one needs to take the 
dynamic, |r2„,.|/(3-parts of the two out of n bubbles, and set = 0, Q — *■ in the remaining n — 2 ones, as any extra 
factor of as well as any extra factors of Q eliminates the logarithmic divergence of the momentum integral. The 

dynamic terms can come either from the two bubbles, already present in the skeleton second-order diagrams, or from 
the bubbles associated with the vertex corrections. It is intuitively plausible that, once the two dynamic bubbles 
are chosen at the n— th order, while the rest of diagram is evaluated at Qm — 0, thus constituting the n— th order 
correction to the static vertex. In other words, it is plausible that the non-analytic, T'^-term in the thermodynamic 
potential can be expressed as the two second-order skeleton diagrams, in which the wavy lines are replaced by exact 
static vertices, r'^(k, — k; k, — k) and r'^(k, — k; — k, k), and bare fermion lines are replaced by solid ones. 

This conjecture, however, needs to be verified, as diagrams for the thermodynamic potential contain combinatoric 
factors. A priori, it is not clear whether these factors, combined with those associated with the selection of the two 
dynamic bubbles, would give just the right combinatoric factors to form the perturbative series for the static vertices. 
In order to verify this conjecture, we explicitly evaluated the T'^-term in the thermodynamic potential to the third 
order in U{Q), and compared the result with that obtained by inserting renormahzed static vertices into two-bubble 
skeleton diagrams. This derivation is presented in Appendix^ We find the two expressions, obtained directly and 
via skeleton diagrams, are identical. We did not attempt to prove that this equivalence holds to all orders in U{Q), 
but the agreement between the two third-order results is a promising sign. In what follows, we assume that this 
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agreement holds to all orders in U{Q), i.e., the non-analytic part of the thermodynamic potential is a product of two 
(renormalized) dynamic bubbles and two exact static vertices, r''(k, — k; k, — k) and r'^(k, — k; — k, k). To first order 
in U, these vertices reduce just to U{0) and U{2kF), respectively. 

Vertices r'^(k, — k; k, — k) and r'^(k, — k; — k, k) are exact in a sense that they include all static corrections, coming 
from the states both away from and near the Fermi surface. The latter produce powers of the static bubble, Umi^m — 
0, Q — > 0) = —Z'^m* /2tt, which, we remind, comes from the states in a narrow range near the Fermi surface In 
other words, the vertices include all corrections, except those coming from the dynamic part of the polarization 
bubble. In conventional notations jJiQ, vertices r(k, — k; k, — k) and r(k, — k; — k, k) are related to r'^(0), defined by 
Eq. H1.16|l . As the incoming momenta are nearly anti-parallel to each other, angle 9 in Eq. H1.16|l can be put equal to 

TT. 

Vertex r'=(7r), as a tensor in the spin space, can be represented as 

r^^p^^siir) ^r\k,-k;k,-k)S^^Sps -T''ik,-k;-k,k)do,sSpy. (4.48) 

Quasi-particle Z-factors, resulting from the self-energy insertions into the fermion lines of the bubbles, can be incor- 
porated into a relation between F*^ and the quasi-particle scattering amplitude 

U0,ys{7T) = Z'rip^^si^). (4.49) 

Next, representing the scattering amplitude in terms of its charge- and spin components 

7T 

/o!7,/35(k, -k) = fa-r./dsiT^) = — [fc{'n')SatSi3S + /s (vr)(TQ-y CT/J^] 

= ^ [{/c(7r) - fs{7T)} So^^SpS + 2MTT)Sa.5Spj] , (4.50) 

and comparing Eq. 14.50() with Eq. 14.48(1 . we obtain 

Z'r'^ik, -fc; fc, -fc) = ^ IUtt) - Un)] , Z^r'^ik, -k- -k, k) = -2^/,(7r). (4.51) 

m* m* 

Substituting Z'^T''{k, -k, k, -k) instead of U{0) and ^■2F(k, -k; -k, k) instead of U{2kF) into Eq. we obtain 

the final form of the non-analytic part of the specific heat in a generic Fermi liquid 

6C{T)/T = --^f^ [f^in) + 3f!in)] T. (4.52) 

On the other hand, 6C{T) does not have a simple closed form in terms of the Landau interaction function, F{9). 
Indeed, the Landau function is related to vertex F'^, defined in Eq. p. 17(1 . rather than to F'^: 

Fa.0,^5i0) - Z^V^^^^sid) = [{7e(7r) - 7,(7r)} So^^Sps + 2j47r)SasSfiy] . (4.53) 

A pcrturbative expansion for F'^(7r) includes static vertex corrections from the states away from the Fermi surface 
but not in its vicinity, i.e., it neglects the vertex corrections associated with n(i7 = 0,(5 ^ 0). Whereas there 
is a simple relation between the partial components of functions T''{9) and F'^(0) 0, in order to relate these two 
functions at a given angle, e.g., = tt, one has to invoke an infinite number of partial components. In this respect, 
the universal sub- leading term in the specific heat is different from the leading, Fermi- liquid term, C{T)/T =const, 
which is expressed in terms of {Fc {0) cos 9) , where Fc is the charge component of the Landau function. 

To emphasize an absence of a simple relation between r''{9) and F'^(0), we present the relation between the charge- 
and spin components of scattering amplitude, /a(7r), and those of the vertex, 7a(7r) (a — c,s), to third order in U{Q). 
Using results from Appendix IfI we obtain 

/e(7r) = 7,(7r) - (7,(7r) - 7,(^)) (7,(^) - 7,(^) + 2(7,(0))) - (7.(0)%(^ " d)); 

fsijr) = jsiTT) - {ls{Ohs{n - 9)), (4.54) 

where (...) stands for averaging over 9. To the same accuracy, the terms of order (F")^ can be expressed via the 
bare interaction potential U{Q = 2fcFsin6'/2) using the first-order relations: F^(6') = f/(0) - U{9)/2 and T'^{9) = 
— {1/2)U{9). The relations between F'^ and U{Q) to second order are presented in Appendix IfI 

The rather complicated relations between Fj (tt) and F^ (tt) are simplified if bare interaction U{Q) is strongly peaked 
at Q = 0, so that F"(7r) is much larger than F^(^^) for a generic 9, including 9 — tt. In this limit, pcrturbative series 
for the relation between Fjf and F" can be summed up exactly, and Fj?(7r) can be expressed in terms of F"(7r) as 



(4.55) 
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In addition, if the condition 73(6*) is met then,, the contribution to the specific heat from 7c (tt) dominates, and the 
singular term in the specific heat becomes 

47r \I + 7c(7r)/ Ep 



3. Non-analytic term in the specific heat for the Coulomb interaction 

The hmiting case of 7c (tt) — > cxd in Eq. (|4.56l) corresponds to the Coulomb interaction in the weak coupling limit, 
when the dimension-less gas parameter, Vs = V2me^ /kp, is small. [We remind that Eq. (|4.5t)|l is valid only in the 
limit 7s(0) <C 1.] [Recall that T'^{0) is expressed via U{0) and U{9), whereas r^(6') is expressed only via U{d) and, 
for generic 0, is proportional to r^.] We see that in this limit r|f (tt) drops out from Eq. H4.56|) . and the singular term 
in the specific heat becomes 

.C(T)/T^-^|-. (4.57) 

For the Coulomb interaction, the T^— term in the specific heat is universal, i.e., it is independent of Ts at small rg. 
This result is in full agreement with Ref. [Blj . 

For the sake of completeness, we explicitly calculate the specific heat for the Coulomb potential by summing-up the 
RPA (ring) sequence of diagrams for the thermodynamic potential (diagrams 1(a), 2(b), 3(c) ... in Fig.!?)). A sum of 
these diagrams results in a familiar expression for S 

2 1-7 (2^)2 ''i-2y(g)n„(f^„,g)' ^- ' 

where V{Q) — 2Tre^ /Q and a factor of two in front of 11^,1 comes from the spin summation. Replacing the Matsubara 
sum by a contour integral and using Eq. (|4.2|) . we obtain the following expression for the correction to the specific 
heat: 

SC(T)/T^- ^ dn ^ / dQQ arg — — jy- -. (4.59) 

The pole of the effective interaction at fl^ = il^ (Q) = VpKQ/2 corresponds to a 2D plasmon mode, where k = 2me^ 
is the screening wave- vector. Near the plasmon pole, 

~ (4^60) 



i-2V(0)n«(n,0) (n + io+)'-n2(<3) 

the argument of the effective interaction changes by — tt when Q intersects the plasmon curve. The plasmon's 
contribution to SC (T) comes from the region 2i}'^/vpK < Q < il/vp. As typical ~ T , one can neglect the 
contribution from the lower limit of Q for T <C KVp. We then obtain a universal and interaction-independent plasmon 
contribution to dC{T) 

Cp,(T)/T.-^^. (4.61) 



For the contribution from the particle-hole region, we find 

dTSn^T^Jo smk^n/2TJn/v 



Cph{T)/T = - — -^ dn—^—— I dQQtan-'W{Q,n), (4.62) 



where 



w{Q,n)^ \ (4.63) 

VVpQ^ - 02 (g + K) 
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The integral over Q in Eq. (|4.62() diverges logarithmically at the upper limit, where tan""'^ W {Q, f7) « W (Q, il). The 
divergence is cut &t Q ~ kp as Eq. H4.63|l is only valid for small Q. Subtracting off and adding W {Q, fl) in Eq. 14.62|l . 
we split CpH (T) into two parts as 

CPH (T) - (T) + (T) ; (4.64a) 



The integral over Q for Cp^j| (T) is convergent at the upper limit, and the integral can be extended to infinity-thus 
this contribution is universal. One can readily make sure that typical Q ~ Q/vp — T/vp <ti k, so that n drops out 
from function W {Q, il) in Eq. (|4.63|l . Evaluating the integral, we obtain 

C^S(r)/r=^^, (4.65) 

which differs by a factor of (—1/2) from the plasmon contribution Eq. 1)4.61(1 . The sum of the two contributions is 

SC (T) /T - CpL (T) /T + (T) /T = _^!!^^^ (4.66) 

47r £!jp 



which coincides with Eq. H4.57|l . 



Finally, the second term in Eq. (|4.64a|l . Cp^ (T) , gives a regular, linear- in-T, correction to C (T) . To logarithmic 



accuracy, 



and 



f dQQW {Q,n) (nn/vp) f dQ/Q = {Kn/vp)lnkF/K, (4.67) 



(T) = -Cfg {T) — r, I In r, I . (4.68) 

TT 

We remind that our treatment is applicable for ^ 1. Notice that the r^llnrsl correction to C (T) /T can be 
interpreted as the result of the mass renormalization 

\/2 

Sm/m — TglnrJ^. (4.69) 

TT 

This last expression coincides with that obtained in Ref. -XM by evaluating the low-energy asymptotic form of the 
self-energy. 



V. CONCLUSIONS 



This paper presents a detailed perturbation theory for interacting fermions in 2D and analyzes non-analytic cor- 
rections to the Fermi-liquid behavior beyond the second-order in interaction. We derived a full expression for the 
fermion self-energy near the mass shell, valid to an infinite order in a weak short-range interaction. Recent study 
pJJ] found, that to second order in U, the imaginary part of the self-energy diverges as ln|A| upon approaching the 
mass shell, where A ~ oj — ek = 0. Following this lead, we demonstrated that beyond second order, divergences 
become of the power-law form in |A|, and get stronger with increasing order of the perturbation theory. We identified 
the divergent contribution as originating from the interaction between fermions and zero-sound collective excitations. 
In the perturbation theory, the collective mode coincides with the upper edge of the particle-hole continuum. This 
degeneracy causes divergences at any finite order of the perturbation theory. We demonstrated that a re-summation of 
the power-law divergent terms to all orders in the interaction eliminates the power-law divergences in the self-energy, 
as the zero-sound mode splits off the continuum. The still remaining logarithmic divergences near the mass shell are 
eliminated by a finite curvature of the fermion spectrum. A fully renormalized self-energy is then completely free from 
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divergences. We found that for a contact interaction, the real part of the self-energy vanishes on the Fermi surface, 
while the imaginary part of S behaves as liP' In Near the mass shell, both ReS''* and ImS^ evolve rapidly as scaling 
functions of variables A/(u^a;) and IS.EpjuP' . The first scaling variable emerges after the re-summation of power-law 
divergent diagrams, whereas the second one describes the effect of a finite curvature of the Fermi surface. We demon- 
strated that the interaction of fermions with the zero-sound mode gives rise to a kink in the spectral function near 
A = This prediction is amenable to a direct check in photoemission measurements or in momentum-conserved 

tunneling between two layers of the 2D electron gases. 

In the second part of the paper, we discussed the non-analytic part of the specific heat: (5C(r) oc T^. We found that 
the collective- mode contribution to the fermion self-energy does not affect the specific heat, i.e., a non-analytic term 
in C (T) is determined only by the perturbative part of the self-energy. This result was also verified by calculating 
the thermodynamic potential directly, in both real- and Matsubara frequency formalisms. We also considered the 
T^-term in the specific heat for a generic Fermi liquid. We showed that it can be expressed in a simple way via the 
spin and charge components of the quasi-particle scattering amplitude at angle = tt between the incoming momenta. 
On the other hand, dC{T) cannot be expressed compactly in terms of the Landau interaction function. Finally, we 
found for the Coulomb interaction, not amenable to a direct perturbative treatment due to infrared singularities, a 
non-analytic T^-term in C(T) is universal and independent of the gas parameter, r^, for small Vs- 
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APPENDIX A: BACKSCATTERING CONTRIBUTION TO THE SELF-ENERGY 



In this Appendix, we present the calculation of a non-analytic part of the self-energy, resulting from backscattering 
processes at T = 0. This part of the self-energy contains two contributions: from processes with small momentum 
transfers ki ~ k[ , k2 ~ k'2 [see Fig.^b)] and from processes with momentum transfers near 2kp: ki « ~k[, ^2 ~ — ^2 
[see Fig.^c)! (52- and 2fcF-processes, correspondingly). The contribution to the self-energy from the g2-process was 
obtained in ,13i] . and here we just cite the result. In Matsubara frequencies, 

„,2 



Sg2 i'-^nn k) — —i 



SttEf 



In- 



1 



^7 {^m + i^k) In 



^■m - ^<^k 



(Al) 



We neglected regular terms of order cj,^„ in l|Aip . Although the second term in Eq. (|A1|) is also formally of order 
for a generic ratio of ujm and e/c, we have to keep it in order to compensate for a superficial divergence in the first 
term at tOrn — ±ie/c in Eq. HA1|I . Converting to real frequencies and taking the imaginary part, we obtain 



ImEf ffc,u;) 



SttEf 



In 



Ef 



P' In 



1 



- - - efe) In 



(A2) 



We note in passing that while the logarithmic factors in Eq. IjAlll comes from the integration over boson momenta 
Q that exceeds LOm/vF, the boson frequency (Om) is smaller than the fermion one (wm); e.g., for positive ujm, the 
non-analytic contribution comes from —ujm < < 0. This implies that a non-analytic part of the self-energy cannot 
be obtained within a renormalization-group scheme, in which internal energies in the self-energy diagram are assumed 
to be larger than the external ones. Notice also that the self-energy in Eq. (|A2I) is regular at w = ±efc, since the 
superficial divergences in the second and third terms cancel each other. Both of these terms are of order uP for 
a generic ratio of uj and e^, and hence, to logarithmic accuracy, one can neglect them compared to the first term. 
Eq. ljA2|l then simplifies to 



ImS«(u;,fc) 



(A3) 



87r£'F 

Next, we calculate the contribution of the 2fci?-scattering to the self-energy. Substituting Eq. H4.27|l for the polarization 
bubble near Q — 2kF into the second-order expression for self-energy 

d^Qdn 



E(fc, w™) 



(27r)= 



-G(w„ 



an,k + Q)n„(r! 



(A4) 
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and expanding the quasi-particle spectrum Ck+q near Q = 2kp as 

Ek+Q -cfe + VEkpe"^ + {2k%/m) q, (A5) 
where q = {Q — 2kp)/2kp and ir— 9 is the angle between k and Q {\d\ ^ tt), we obtain, after re-scahng the variables, 



T,2kFii^rn, k) = - 



mU kf 
27r4 



dq 



q^ 



1/2 



dO 



where LUm = ^^m/'^EF , — i^m/^Ep, and ek = ^k/^Ep- Using 

dz vr sgnRe(a) 



+ a 



we can re- write Eq. 

S2fei,(Wm, fc) 



U as 

4-^-3 



1/2 



sgnRe(A/ q- ek - i{uj,n + ^m)) 



(A6) 



(A7) 



(A8) 



As we are interested in a Inw contribution, we need to check the behavior of the integrand for large |g| (logarithms 
come from the integration over the range where internal variables are larger that the external ones). For q — > +00, 
the integrand approaches a finite limit (\/2) as \/2(l + ^^m/^?^); neither the leading nor sub-leading terms in this 
asymptotic form produce a logarithm upon integrating over q. However, for negative and large g, (q-t- \J q^ -\- H^)-'"/^ « 

\q\ sgn(n„ 



and \J q — tk — i(i^m + ^m) ~ — ^-v/M" sgn(r2m +Wm), so that the integrand behaves as \q\ ; hence 
the logarithmic singularity does come from this region of q. To logarithmic accuracy, we obtain 



Evaluating the integrals, we find 



4-^-3 



sgn(f2^ 



-Ef/vf 1 



T,2kF{^m,k) = -i— In 

OTT hp lUJ. 



Ep 



Continuing to real frequencies and taking the imaginary part yields 



(a;,fc) 



Ep 



" 2 1 
iu In 

SttEp \uj 



(A9) 



(AlO) 



(All) 



Comparing this result with Eq. ljA3|l . we see that the non-analytic parts of ImS^ resulting from the g2— and 
2fcj?— processes are identical within logarithmic accuracy. The result for the sum of the two terms is quoted in 
Eq. 1)3. 5|l . Observe that, according to Eq. HAlip . a non-analytic, w^ln|w| part of T,2kF{^m,k) comes from negative 
q, i.e., from Q < 2kp, where the static polarization bubble is just a constant; thus a singular correction is entirely 
dynamical. A singular, {Q — 2/?^)^/^, behavior of the static bubble n(f7 = 0,(5) for Q > 2kp does not give rise to 
an imaginary part of the self-energy- this result follows due to the fact that a static density fluctuation cannot decay 
the quasi-particles. However, the fact that even for a dynamic bubble, only the region of Q < 2kp is responsible 
for a non- analytic part of S(fc,w) is rather peculiar and has not been emphasized explicitly in earlier work jl3l |. 
Note in this regard that both regions of Q < 2kp and Q > 2kp contribute to a non-analytic behavior of the spin 
susceptibility 0, ^ . 

Now, we evaluate the integrals in Eq. I|A8|I beyond logarithmic accuracy. The result is 



ImS^,.^(fc,o;) 



8nEp 



2 1 

J In - — - 



4n 



- {uj + Cfe) In 



fk 



£k 



(A12) 



Again, the last two terms are of order uj'^ for a generic ratio of to and but we keep them in order to demonstrate 
that ImE^j.^(aj, k) remains finite at uj — ±6^. 
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APPENDIX B: NON-SYMMETRIZED VERTEX 



In this Appendix, we derive an expression for forward-scattering part of the non-symmetrized vertex, F, summing 
up diagrams with the maximum number of polarization bubbles to all orders in contact interaction, U. We then anti- 
symmetrize the vertex, and substitute the result into the Dyson equation [ Eq. (|3.2U|) ] to obtain the corresponding 
part of the self-energy. The diagrams for a non-symmetrized vertex up to the third order are presented in Fig. |3| In 
the Matsubara technique, we associate a factor of —U with each of the interaction lines, and a factor —2 with each of 
the polarization bubble. There is also an extra factor of —1 for exchange processes in which the two outgoing legs are 
permuted (the last diagrams of second and third order in Fig.|31l. We present here a general recipe for calculating the 
j^th y j^-j Qj.(jgj. vertex diagram. As is seen from the Fig.|2| the vertex consists of two parts. The first part comes 
from the direct interaction and contains a spin factor Sa-yS/^s- The second part is due to the exchange interaction, 
and comes with a spin factor SaeSp-y- At each order, there is only one exchange diagram whose contribution is 
— {—UYll''~^6aeSi3-f. (At second and third orders, these are the first and third diagrams in the third column of Fig.O 
respectively.) The rest of diagrams are due to the direct interaction and contain various number of bubbles. At order 
I/, the number of these bubbles (i?) varies from to i/ — 1. For a diagram with R bubbles, R + I interaction lines are 
used up in making a chain of bubbles and connecting it to two external solid lines. The remaining N = v — R — 1 
interaction lines can be arranged anywhere either at the two ends of the chain of bubbles or inside the R bubbles. 
There are 5* = i? -I- 2 sites where N interaction lines can be placed. The number of diagrams with R bubbles is equal 
to the number of ways to arrange N lines among S sites: 



{S + N -1)1 



(s'-i)!iv! {R + ly.iiy - R - ly: 

Consequently, the contribution to F from diagrams with R bubbles is 



(Bl) 



{-Uy{-2fn'^-^S^y6f3e. (B2) 



{R + mu-R-iy.' ' ^' ""^"'^^ 

The total contribution from all bubble diagrams at the order v is then 



E 



{R + mv-R-iy: 



{-UYTl'^-'5^^5p, = - ^ \^^u''W-'5o.y5p, (B3) 



where we have used an identity 



E 



= (i-2r = (-ir. (B4) 



Adding up direct and exchange terms, we obtain the following form for non-symmetrized vertex at order v > 1: 

Kpae = -^(1 - {-rW^^^'S^ySp, - {-UyW-^S^.Sp^. (B5) 
The vertex function can now be readily summed up to all orders, with the result 

u , , c^'n(g) 

Using an SU(2) identity 

5o.e5p^ = (1/2) {<^'y% + 5o.y5p,) , (B7) 

and introducing dimensionless spin and charge vertices Qp and spin Qcy, defined in Eq. (|3.16ll . we obtain Eq. 13.15|1 . 
Anti-symmetrizing the vertex as prescribed by Eq. (|3.17() and substituting the result into the Dyson equation 13.20|l . 
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we obtain for the self-energy 



Sa/3(p) = s^p / uG{p-q)- / c/r^„;^0G(p')n(p-p')+ / UT^^.p^G{p")n{p-p") 

J q J p' J p" 



{irr|nF + r^}<=<''"W''-''") 



V« / VG{p -q]+ V« / "'^^, G(P - 9) + 1 ^^2!l^G(p - ,). (B8) 



Re-arranging the resuh, we obtain the self-energy in the form of Eq. H3.22|l . 



APPENDIX C: RELATION BETWEEN EQS.lHTst AND (IHTll . 

In this Appendix, we discuss the relation between the two results for the specific heat: the one derived in Ref. ^ 
[Ea. H4.5|l ] and the one following from the Luttinger-Ward functional for the thermodynamic potential [Eg. H4. 11(1 ], 
We show that the two expressions yield identical results for the T^-term in the specific heat, provided that one uses 
in temperature-dependent self-energy in Eq. 14.5|l (the recipe in Ref. l] is to use the self-energy at T = 0) . 

For the sake of simplicity, we focus on the contact interaction and consider only a perturbative part of the thermo- 
dynamic potential, neglecting the interaction with the zero-sound mode. 

We begin with the Luttinger-Ward approach. To second order in u, the Luttinger-Ward formula for the entropy is 



5 - 5o + ^ 



Me/eS(w„i, fc,T)G'o(Wm, fc) 



(CI) 



Here 5*0 is the entropy of an ideal gas, A'o = m/2i: is the density of states per spin projection, Gq = {iuJm, — ^k)~^ is the 
bare Green's function, and E(wm, fc,T) is the self-energy, which depends on T both implicitly, via = 7rT(2m-|- 1), 
and also explicitly. The second-order self-energy in Matsubara frequencies is, to logarithmic accuracy, 

E(a;m,/c,r) = -——u^T'S^ |ri,„|sign(w,„ + r2m) In — — • r-TT^w^^V' |rim|sign(wm -|-ri,„) In . 

2hF efc - z(2i2„i + cj,„) 4i/_F ^ ek - i(^m 

(C2) 

The first term in Ea. ljC2|) comes from backscattering and includes the sum of contributions from the g2~ and 2kp- 
processes. The second term in Ea. (jC2ll comes from forward scattering. 

Had there been no dependence on under the logarithm in HC2|I . the integration over in (|C1|I would have been 
straightforward, as it would have involved only Go(tt'm, k). Using the familiar relations 

V J defeReG^(c^, e^) = 0, ^ defeImG«(c^, e,.) = -ir; (C3) 

and converting the Matsubara sums into the contour integrals, one could readily verify that the expression for the 
entropy would have involved only Rel]^(w, fc, T) on the mass shell and that both terms in (jC2p would have contributed 
to the entropy. 

Because of the fc— dependence under the logarithms in HC2|I . the actual situation is different. Indeed, substituting 
the second term in HC2|I into (jCip we find that the integral over efe reduces to 

defe , Ep 1 , 2 Ep — iuj 



In = - In" — ^ (C4) 

-Ef '^k - Efc - «t^m 2 -Ep - lUJ 

The result in vanishes in the limit of \ujm\ ^ Ep. Thus, to logarithmic accuracy, forward scattering does not 
contribute to the entropy. It is instructive, however, to see how the zero for the integral in ljC4|l comes about. In 
fact, there are two contributions to this integral, each of order In Ep /\ujm\- The first comes from the region of small 
ek- e/c — '^m- For this contribution, the dependence on ek under the logarithm can be neglected, i.e., the logarithmic 
factor can be approximated by hiEp/\u!m\- The integration over ej, then yields 

Ep f dek . Ep 
in- r / : = iTTSgnw™ In -j -. (C5) 
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The second contribution comes from the region of large Cfe : efe ^ Wm- Here, we obtain, combining the contributions 
from positive and negative 



lUJr, 



= -i7rsgnw,„ln-— ^. (C6) 

Combining IIC5I) and l|C6|l . we find that for the forward scattering contribution to S, the net integral over vanishes. 
For the backscattering contribution, the integral over is 

In , J' ... . (C7) 



There are again two contributions to this integral to order InEp /\(^m\- The contribution from efe ~ w„i is the same 
as (|C5|) . while the one from Cfc is 

— iTTSgnfJ,,, In -j r . (C8) 

It follows from ljC2|) that, at low temperatures, the dominant contribution to the self-energy comes from such Vim 
that are of different sign compared to Wm- As a result, for the backscattering part of the self-energy, the small-efc— 
and large-e^-contributions to the integral over add up instead of canceling out. 

The outcome of this analysis is that, to logarithmic accuracy, one can still formally neglect the dependence on 
Efe in the logarithm in (|C2p : however, one should simultaneously neglect the forward-scattering contribution to the 
self-energy, and multiply backscattering contribution by a factor of two. Converting the Matsubara sum in Eq. (jCl|) 
into a contour integral, one can express the entropy via the real part of the self-energy on the mass shell, as 



S^Sq-No du 



(C9) 



We emphasize again that this expression is two times larger than the one that one would obtain by just neglecting 
the dependence on Cfe in the forward scattering self-energy. 

Next, we find the relation between the two terms in (|(y9(l . If the self-energy were independent of T, the second 
term would not contribute, and the entropy would be given just by the first term. We show, however, that the two 
terms in ljC9p are in fact equal, i.e., the entropy can be formally re-expressed via only the first term in ljC9|l . but with 
another extra factor of two. 

To demonstrate this, we need an expression for ReS^(ci;,r). The easiest way to obtain it is to convert Matsubara 
self-energy into \mTi^{Lo,T) and then obtain Rel]^(ci;,r) by a Kramers-Kronig transformation. Approximating the 
logarithm in the first term in ljC2p as XiiEp /\Vlm\: we obtain 



,,2 foo 

Iml]^(w,r) = / dfifiln 

-iTTEp J-oo 



Ep 

m 



, n , oj + n 

coth tanh 

2T 2T 



(CIO) 



Substituting the Kramers-Kronig transform of IjClOp into ljC9|l , and integrating in the second term in ljC9|l by parts, 
we find after some algebra that the two terms in ljC9p contribute equally to the entropy. 
As a result, the net expression for the entropy in the Luttinger-Ward formalism is given by 



oo 



S = So-2Nol du;^^ei:"iu;,T). (CU) 



oo 



It differs by a factor of four from what one would have obtained neglecting both under the logarithm in ljC2|) and 
the T-dependence of ReE^(w, T). 

Eq. (|(T1|I is very similar to the one from Ref. which, to the lowest order in the interaction, reduces to 

S^So + ^N„ p duj I dek^^lm[j:{Lu,k)Go{u,ek)]. (C12) 

We first assume, and then verify, that once e^.) is evaluated in real frequencies and substituted into (|C12|I . the 
momentum dependence of the self-energy can be neglected. Using (jC3p . we then obtain that (|C11|) and ljC12p coincide 
if one substitutes the T -dependent self-energy in ljC12|) . 
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Finally, we show that the momentum-dependent part of S(a;, k) does not contribute to the entropy. The proof is 
somewhat tricky. The backscattering term in the second-order self-energy [first term in HA2j does indeed contain the 
dependence on under the logarithm, just like in the Matsubara self-energy. Because of this dependence, there is 
a non-zero contribution to the entropy from ImE-'^ReG^ term in (|C12p . Should one keep it, this would modify the 
answer for 5 by a factor of two, compared to (|C11|I . 

However, in contrast to the Matsubara self-energy, which is strictly perturbative, the self-energy in real frequencies 
does develop mass-shell singularities in the forward-scattering part. We remind that the reason for these singularities 
is the presence of the zero-sound mode that, within the perturbation theory, coincides with the upper boundary of 
the particle-hole continuum. These mass-shell singularities give rise a non-perturbative part of the self-energy, Szs 
and also modify the part of the self-energy coming from the interaction with the particle-hole continuum, Sox + 
[see Eq. (12231. 

We show that the two contributions to the entropy - one from ImS^ReG^ and the other from ImS^ReG^ - cancel 
each other. To this end, we demonstrate that 

2t: J w - efc 

where Y,^ is a sum of backscattering and forward-scattering contributions, does not contain terms non-analytic in ui . 

We begin with the backscattering part. As it was discussed in Sec. there are two types of backscattering 
processes: 172 and 2fci?. The self-energy from the g2-pi'ocess is given by Eq. ljA2|l . Terms that are independent of ek 
obviously do not contribute to B (uj) . Substituting the rest of Eq. ljA2|l into (|C13|I . we obtain 

2 

m 11 

S,.=-^;^c.[c.[. (C14) 

In (jC14p . we neglected terms of order lu, which contribute just to the renormalization of the effective mass at T = 0. 
It can be also verified that the subleading term to (|C14p is analytic and scales as uj^/Ep- 

The 2fci?-contribution to the self-energy is given by ljA12ll . Substituting this contribution into ljC13l) and evaluating 
the integral over ek, we find that i32fef.(w) contains regular terms of order w, oj^, etc., but no w[ct>[ term. Therefore, 
for our purposes, 

B2kA^)=0. (C15) 

Consider now the forward-scattering (174) part. The imaginary part of Ezs is analytic near the mass shell and does 
not contribute to B. The imaginary parts of SpH and Sox, however, depend on A = w — Cfe logarithmically away 
from the immediate vicinity of the mass shell, where the logarithmic dependence is regularized [see (|3.26|l . (|3.31|) and 
(|3.49|) ]. Substituting the regularized expression, Eq. H3.49|l into HC13|) we find that the Bg^ reduces to 

^-=16^^ H-y_;^ln|.|^^^^|c.|c.y^ (C16) 
where (5 oc w is positive. Performing the integration, we obtain 

777 II 

B9.^^^^H, (C17) 

We see that it is opposite in sign to the contribution from the 52— process, Eq. (|C14|I . Adding up the three contribu- 
tions, Eqs. ljC14|) . I|C15|) and (|C17|I . we obtain that the total B{lo) = 0. This proves that the momentum dependence 
of S in HC12|I can indeed be neglected. 



APPENDIX D: SPECIFIC HEAT IN ELIASHBERG-TYPE THEORIES 



In this Appendix, we consider the specific heat for the case when the fermion self-energy depends only on the 
frequency but not on momentum. Such a situation occurs, e.g., in Eliashberg-type theories, which describes the 
interaction of fermions with slow boson modes We show that in order to obtain a correct form of the linear- in-T, 
Fermi-liquid contribution to the specific heat, one can use the approximate relation Eq. (|4.1|l (from Ref. Q), which 
expresses the thermodynamic potential solely in terms of an exact fermion Green's function. However, a correct 
form of the sub-leading, non-analytic part of the specific heat can only be obtained using the full Luttinger-Ward 
expression, while Eq. H4.1|l gives an erroneous result, even if one uses the self-energy at finite T, Yi{ujm, T). 
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A general form of the thermodynamic potential, S, for a system of interacting fermions is given by Luttinger-Ward 
formula, Eq. H4.8|l . It expresses S in terms of the full Green's function and infinite series of skeleton self-energies. 
In many cases, though, the interaction between fermions is strongly enhanced in a particular interaction channel. In 
such a case, multiple interactions between fermions in the same channel can be adequately described as an exchange 
of corresponding low-energy bosonic collective modes. If, in addition, these modes are slow compared to fermions, 
the self-energy of fermions will be independent of the momentum, I](cjm,fc) = and corrections to fermion- 

boson vertex can be neglected due to the Migdal theorem. For these cases, the Luttinger-Ward formula reduces to a 
closed-form expression in terms of full propagators of fermions and low-energy bosons 



In [el + {Um + - iT,{uJra)G{k, COm) 



[HD-\Q, + 2n(r!„, Q)D{^,n: Q)] 



(Dl) 



the full fermion Green's function, _D(Q,il„j) 



Here G{k,Um) = [i {w„ + i;(w„J} - efc] 

{D^^iQ, n) - 2n{n,n, Q)y is the fuU boson propagator (a factor of 2 is due to the spin summation), Dq{Q,VI) 
is the bare boson propagator, S(f2„j) and n(ri,„,(5) are the fermion and boson self-energies, respectively, and g is 
the effective fermion-boson coupling. [In this Appendix, the Matsubara self-energy is defined with an i up front.] 
Eq. jnil) was first obtained in the context of the clcctron-phonon interaction 0, l44l | (in which case is the phonon 
propagator), and was applied later to the electron-electron interaction mediated by Landau-damped collective 
modes ^3 Ulii 4L 4S]- In the latter case, the frequency dependence of D{Q,Q,m) comes predominantly 
from the boson self-energy, and the bare boson propagator D(){Q,nm) can be approximated by its static form 
Dq (Q) = Do{Q,0) 0,ESE3iE3- Notice that Eq. I|D1(I is applicable to the interaction with both charge- and spin 
modes, except for the spin case it has to be modified slightly due to the spin structure of the interaction |49l |. 

The fermion and boson self-energies, E(a;m) and n(r2,„), arc obtained from the condition that S is stationary with 
respect to variations in E(wm) and n(f2„i). Conditions 



SE/SY.{k) ^ SE/Sn{k) = 



yield |44 



j:M = tTg'Yl J G(fc',c.;„)i^(k-k',c 
U{Q^,Q)^-Tg^J2 I S^G{k,u^)G{\<i-Ql,Um-^m)- 



Using Eq. (|D2|I . the last term in Eq. ljDl|l can be re- written as 



'^E 



(2^) 



D 



n(r!„)7^(Q, r! 



or, equivalently, 

Accordingly, Eq. HD1|I reduces to 

S = -2T5: 



2^^E / (0j5](c.„)G(fc,..„). 



(D2) 



(D3) 



(D4) 



(Pk 



2 



(2^)^ 



k + {^m + S(w,„)}^ - iE(cj,„)G'(fc, UJra) 



(D5) 



43 



or, equivalently, to 



In 



[\n[D-\Q,n^)\^2I].{n^)D{Q,^r. 



(D6) 



Each of the last two expressions can be simphfied further. In Eq. IjDSp . we switch from the integration over momentum 
to that over usin g / d^k/{2n)^ — Nq J dek, where A^o is the density of states. Integrating over e^, we find that 
the first Une in Eq. I|D5|I reduces to 



(D7) 



As the sign of S(a>„i) coincides with that of cj,„, the self-energy drops out, and Eq. I|D7|I reduces to the result for 
non-interacting fermions, SpQ. Substituting Eq. (|D7p back into Eq. (|D5(I . we obtain 



where 



d^'Q 
{2t:)D 



ln[D-\Q,n„,)], 



Efg^-TNoJ2\ 



(D8) 



(D9) 



is the thermodynamic potential of a free Fermi gas. Eq. (jPSp is often called a Luttinger-Ward expression for the 
thermodynamic potential in Eliashberg-type theories. For a spin interaction, the factor 1/2 in the second term in 
Eq. (|D8|1 is replaced by 3/2 (491]. We note in passing that the frequency sum in Spc formally diverges, but its 
temperature-dependent part, which is what we need, can be extracted by using the following spectral representation 



dx 



The above relation in conjuction with the following identity 



1 1 f . 



(DIO) 



(Dll) 



can be used to extract the T-dependent part of Spc: we thus obtain the familiar Fermi-gas result, Cfg — (27r^ A*'o/3)r. 

For the thermodynamic potential in the form of Eq. (|D6|I , we can use the fact that for Landau-damped collective 
modes 11(0^) c« \^m\ and expand in 11, as higher powers of Vl in the summand in Eq. (|D6|) generally lead to higher 
powers of T in S. Expanding the logarithm, we see that the term linear in 11 drops out, and the thermodynamic 
potential is given by 



(fc,c.) 



(D12) 



This is the same expression as Eq. (|4.8ll , obtained by a different approach as compared to that in Ref. 1] . In a Fermi 
liquid, S(cL'm) = XiOm at the lowest frequencies. Substituting this form into Eq. ljD12|) . we immediately obtain 



Cfl = Cfg(1 + A), 



(D13) 



Thus, both approaches-the one based on the Luttinger-Ward functional and the one used in Ref. Q-give the same 
result for the FL-part of C{T). 
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The issue now is whether the approximate form of S [Eq. (|D12|l ]. when modified to include a finite-temperature 
self-energy can correctly describe the non-analytic corrections to the Fermi liquid. We argue that it does not. Indeed, 
expanding the logarithm to second order in Hm in Eq. l)D6p . we obtain 



-T 



(DM) 



To avoid further complications with a long-range interaction, we assume that Dq{0) is finite. Then, as we showed in 
Sec. lIVBl the integral over Q gives In \ ^lm\ which, upon summation over f2m, results in a non-analytic term in S 
{T^ in 2D and T^lnT in 3D-see Appendix G). For these situations, the boson contribution to Eq. (|D6p cannot be 
neglected when evaluating the non- analytic term in the specific heat. 



(D15) 



We note in passing that expanding the logarithm in Eq. ljD8() to order H Dq results in 



" = ^FG 



,n\q,n)Dl{q). 



Comparing Eqs. ljD14|) and IjDlSj) . we see that the first term in Eq. (|D14p must be twice the integral of IPDq. This is 
consistent with the observation we made in Sec. lIV Al where we found a factor of four difference between — 2TlnG'~^ 
[the first term in Eq. (|D14|) ] and Eq. I|D15|) . However, a factor of two difference was attributed to neglecting the 
temperature dependence of S](a;,„) while converting from Matsubara to real frequencies. 

For completeness, we also demonstrate how one can obtain the Fermi-liquid result, Eq. I|D13|I from Eq. I|D8I) . which 



expresses in terms of the boson propagator. As we have already mentioned, in order to get a Fermi-liquid, T -term 



in S, one has to expand to first order in n(ri,„, Q) cx \ flm\- To reproduce Eq. (jD13|) . one therefore needs to relate the 
Landau damping term to A. This relation can be found for arbitrary Dq{Q). To shorten the presentation, we just 
consider a model form of Dq {Q)-a, Lorentzian peaked at Q = 



DoiQ) 



Do 



For this form of Dq{Q), the fermion self-energy at the lowest frequencies is readily obtained from Eq. 1D2|I : 



(D16) 



The polarization bubble at low frequencies and small momenta is also obtained from Eq. (|D2 



n(r!™,Q) = 



7 l^ml 

^0 Q ' 



(D17) 



(D18) 



where 7 = mg'^ / (nvp) (we set interatomic distance a = 1). Substituting Eq. IID18|) into Eq. I|D8|) . integrating over 
boson momentum and collecting terms of order in S, we obtain 



(D19) 



Summation over boson frequencies is performed by using spectral representation, in the same way as the sum over 
fermion frequencies in Eq. ljD9(l . For temperature-dependent parts of the two sums in Eq. I|D19|I we find 



Substituting this into Eq. (jD19|) . we obtain 



rj-iZ 



Comparing Eq. (jD2ip and Eq. HD17|I , we see that 



1 + 



9 



(D20) 



(D21) 



[1 + A] 



(D22) 
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i.e., Eq. (jPlSp is reproduced. 

Finally, we discuss the specific heat near a Quantum Critical Point (QCP) which, formally, corresponds to the limit 
of ^ = oo in Eq. (|D16|I . Here, we find 

S = Spo + ^7^/^T^|f7™p/^ (D23) 

Evaluating the sum using 

where u = dz/{z^ + 1) = 27r/(3%/3), we obtain 



S = Sfg - ^^^y^T'/'. (D25) 



Using Eq. (|4.2|l . we obtain 



C(T) = Cfg + ^^1^72/3^2/3^ 

TT 

Eq. (|D25|I also allows one to verify the conjecture in Ref. 0| that Eq. ljD12p can be used to evaluate the specific 
heat beyond the Fermi-liquid term. Evaluating the self-energy at QCP, we find jT^. li^ 

•^f \ 2/3 1/3 1/3 1 /T^O'7^ 

27rV3 fF7 ' 



Substituting this result into Eq. (|D12|I . and evaluating the momentum integral and the frequency sum, we obtain 

i^fc , r 9 . on _ 0.3546 



(D28) 



which differs from Eq. (|D25p by a numerical prefactor. We see that using the zero-temperature self-energy in Eq. HD12p . 
we do not reproduce Eq. I|D25(I . This is another indication that Eq. l|D12p is not valid for the calculations of C(T) 
beyond the Fermi-liquid term. 

To avoid confusion, we emphasize that Eq. ljD6|) is valid only for Eliashberg-type theories. For the problem that 
we considered in the main text, the self-energy is /c— dependent, and this eventually makes Eq. 14.8|1 valid to second 
order in U, provided that one uses a finite-T self-energy instead of a zero-T one. In view of the above consideration, 
however, we do not expect Eq. H4.8|l to remain valid at higher orders in U. In any event, it is always safe to use 
Eq. IjDSjl in the calculations of the specific heat. 



APPENDIX E: EVALUATION OF THE THERMODYNAMIC POTENTIAL IN REAL FREQUENCIES 



In this Appendix, we find the specific heat by calculating the thermodynamic potential which is expressed in real 
frequencies. The evaluation of the second-order contribution to S is straightforward: we just replace the Matsubara 
sum in Eq. (|4.33l) by a contour integral. Using Eq. 14.21) . we obtain 



SC (T) /T : 



9 d 

dT 



dn 



47r T2sinh^f7/2TJ {2n 



(El) 



where 11^(51, Q) =lmlP (ft + iO'^, Q) [see Eq. H3.1()|l ]. When differentiating S with respect to T, we assumed that the 
non-analytic part of the particle-hole bubble does not depend on temperature-keeping this source of the temperature 
dependence would result only in analytic, T^-terms in C (T) /T. Using Eq. H3.11|l for the singular part of 112^(17,(5) 
near (5 = 0, performing elementary integrations, and multiplying the result by 2 to account for the contribution from 
2kp, we indeed reproduce Eq. (|4.23l) . 

Next, we demonstrate how non-pcrturbative contributions to C (T) cancel out in the real- frequency formalism. To 
this end, we use series Eq. (|4.25|l for the thermodynamic potential 



-2UU„ 



f-^lnt/p-|ln(-g.) 



(E2) 
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FIG. 9: Argument of charge vertex, (j> = axg Qp, where Qp is defined in Eq. I)8.16|l . as a function of the boson momentum Q at 
fixed frequency Q. Q — Q/c corresponds to the position of zero-sound pole, whereas Q = Q/vf corresponds to the position of 
particle-hole contiimum boundary. 



where the effective vertices of charge and spin channels, Gp and Qa, are defined in Eq. H3.16|l . We recall that the 
charge term -(1/2) In^p -contains the contribution from the zero-sound mode. It is convenient to single out first- and 
second-order contributions in U, i.e., to re-arrange Eq. ljE2p as 



^ — ^0 



^1 i- ^2 -I- ^ 



(E3) 



where 



un, 



(c/n„ 



(E4) 



and 



2S2 



1 



(E5) 



contains a combined contribution of all higher orders. Our goal is to show that there is no non-analytic contribution 
to the specific heat from S'. Converting Matsubara sums into contour integrals, we obtain for the non-perturbative 
part of the specific heat corresponding to S' 



2 arg^^p 



27rT2sinh^r2/2r7 (27r)^ 



-arg(-g,)-C/2nf 



(E6) 



When differentiating Eq. (|E6|I . we assumed again that the only T~ dependence comes from the Bose distribution 
function. The T— dependence of vertices becomes important either near a finite-temperature critical point 8] or in 
ID; neither of the cases are considered in this paper. 

The momentum dependence of aigGp at fixed > is shown in Fig.lHI (solid line). For Q < fi/c, where c is the 
zero-sound velocity, ReQp > and ImQp = 0, hence argC/p = 0. At the zero-sound pole, Q = fl/c, aigGp jumps from 
zero to — TT. It remains equal to — tt until Q reaches the boundary of the particle-hole continuum, Q = fl/vp. At this 
point, arg^p jumps from — tt to — 7r/2 and decreases monotonically upon further increase in Q. On the other hand, 
argC/cr is nonzero only inside the particle-hole continuum {Q > VL/vp). Finally, as oc 5{VP' — VpQ'^), this term is 
only relevant at the boundary of the particle-hole region. This behavior of argC/p, aigGa, and 11^ suggests that it is 
convenient to split the momentum integral into the one over Q < fl/vp and another one over Q > ^/vp, and consider 
three terms in Eq. HE6|I separately. In the region Q < fl/vp, only the collective mode contributes to the momentum 
integral in Eq. (|E6|I . and this contribution is given by 



(fQ 1 



Q<n/vp (Ztt) ^ 



1 

2^ 



n/c 



dQQ 



1 \ n'^ 



(E7) 
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At the last step, we have expanded the full result to order li} . The contribution from arg^p in the particle-hole region 
Q > V,/vp is 



/ 

Jc 



Q>n/vF (Stt)^ 2 



-1 



— tan 



n 



n 



1 + uy 16z;^ 



164' 



(E8) 



To ensure the convergence of the Q-integral at large momenta, we subtracted off a term proportional to Imll^, which 
gives no T^— contribution to C{T), and used that dx (ta,n~^ x~^^'^ — x~^^^) = — 7r/2. The contribution from 
arg Qa- comes only from the particle- hole region and is equal to 



jfQ_3 
>n/vp {2Trf 2 



tan 



n 



^l-U ^4Q2 _Q2j 1 - U ^4Q2 _ Q2 



(E9) 



Again, a term proportional to Imll^ has been subtracted off to ensure convergence. Finally, the t/^n|^-term in 
Eq. ljE6|l . coming from the boundary of the particle-hole region, yields 



1 of?^ 

4 vl 



(ElO) 



Adding up Eqs. ^E7IE10|) . we find that 



C (T) = 0, 



(Ell) 



as it was anticipated. Once again, this means that non-perturbative corrections do not change the result at order 
M^-it is stiU given by Eq. (|4.14|l . 



APPENDIX F: EVALUATION OF THE DIAGRAMS FOR THE THERMODYNAMIC POTENTIAL 



1. second-order in the interaction 



We have shown in the main text that, to second order in U{Q), the thermodynamic potential, S, contains a 
non-analytic, T^-term whose magnitude depends only on U{0) and U{2kF)- To this order in the interaction, S 
consists of two particle-hole bubbles, and the argument for the non-analyticity in S was that it originated from both 
Q = and Q — 2kp non-analyticities of the bubbles. This argument, however, does not specify the relation between 
fermion momenta in the two bubbles, and therefore it does not distinguish between the cases when the total incoming 
momentum in the two vertices of diagrams 2(b) and 3(a) in Fig.0is near zero or near 2kp. 

Now, we look into diagrams 2(b) and 3(a) in Fig. 13 in more detail, and show that the non-analytic term in 
S involves only vertices with "ID" momentum structure (k, — k; k, — k) {Q — contribution) and (k, — k; — k, k) 
{Q = 2kF contribution). For definiteness, we consider diagram 2(b) and focus on the (5 = contribution. The 
2fci?-contribution to diagram 2(b) and diagram 2(a) can be treated in a similar manner. 

The argument why the momenta in the two bubbles are related to each other is based on the observation that in 
order to have S = 5n oc , the summand Sq oc J dP'QlP(Q,m, Q) must be non-analytic in Vl. The momentum 

integral does, indeed, diverge logarithmically at vpQ ^ f^m, as n^(ri,„,(5) contains a term ri^/Q^, which is just 
the product of the flm/Q terms in each of the bubbles. Then Sji oc ri^,j In \ flm\i and the summation over Matsubara 
frequencies yields S oc T'^. 

We now demonstrate that the |r2m|/Q term in n(ri,„,(5) comes from integration over internal momenta k in a 
narrow range around k • Q = 0. To see this, we recall that the bubble has the following form - 

n(f^™,Q) = (l/27r2)T^ / d^kGoiLOm,'k)GoiLu + n,k + Q). (Fl) 
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k + q 




-k + q + q 




k+q+q2 



FIG. 10: A second-order ring diagram for the thermodynamic potential. Momenta qi and q2 are smaU: qi,q2 <S 
Labelhng on the left and right diagrams corresponds to backscattering and forward scattering, correspondingly. 



Expanding Ck+q near the Fermi surface as Ck + vfQcos9 and replacing the integration over (Pk by that over de^dO, 
we obtain 



n(an,0) j deu J dO 



1 



{iuj,-a - Ck) {iuJra + i^m " Cfe " VpQ COS 6) ' 



(F2) 



For the dynamic part of 11 (0^, Q), the order in which the integration over efc and frequency summation are performed 
does not matter. Integrating over first, and then summing over frequency, we obtain 



Integration over 9 gives 



n(r2,Q) (X j dO 



n(an,Q) cx 



i^m — VpQ COS 9 



(F3) 



(F4) 



It is important that typical | cos 6*1 in the angular integral are of order |f2m|/ui?Q. For large Q, Eq. (|F4|) reduces 
to n(rim,Q) oc \flm\/Q, whereas typical angles are near ±7r/2: \9 ± 7r/2| ~ \ilrn\/vFQ, i.e., k and Q are nearly 
orthogonal. 

Since typical momenta in both bubbles in diagram 2(b) of Fig. [7|are nearly orthogonal to the same vector, Q, 
they are either almost parallel or antiparallcl. In the first case, the momentum structure of both vertices in diagram 
2(b) is (k, k; k, k), in the second, it is (k, — k;k, — k). We now demonstrate that only the backscattering vertex, 
(k, — k;k, — k), contributes to the non-analyticity in S. To see this, we evaluate diagram 2(b) in a different way, 
namely, via integrating the product of four Green's functions over a common fermion momentum (momentum k in 
Fig. I1(J|) rather than pairing them into bubbles. As shown in Fig. ^1 the momenta of the four Green's functions 
involved are labelled as either GkCk+qiG-k+qaC-k+qi+qa (backscattering) or GkCk+qjCk+qaGk+qi+qa (forward 
scattering). The arguments presented above for a single bubble suggest that the non-analyticity in S comes from the 
region of small qi and g2- We show that the integral over small qi and q2 indeed gives the correct result, but only in 
the first case, when two fermion momenta are near k and another two are near — k, while in the second case, when 
all four momenta are near k, the integral over small qi and q2 vanishes. 

We begin with the first case. We first assume, and then verify, that one can expand the quasi-particle spectra to 
second order in qi and q2. For a circular Fermi surface, we have 

a : 1i a : 12 fa, a \ + ^2)^ 
Ek+qi = efc-|-u_F(7iCOs6'i + — ;e_k+q2 = ffe -''^f'?2 cos 6*2 + — ; e_k+qi+q2 = efc-w_F((?i cos6'i-|-g2 cos6'2jH — . 

""(FS) 

We also assume, following the reasoning for a single bubble, that both qi and q2 are nearly orthogonal to k. This 
implies that qi and q2 are either nearly parallel or nearly antiparallcl. Collecting the overall factor for diagram 2(b), 
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expanding in 9i and 62 near ±7r/2 and introducing new variables x = vpqiCosOi, y = vpq2COs92, we obtain the 
contribution to diagram 2(b) from small qi and (72 



j-2 ^Q'^ COQ poo poo poo poo 

(27r)5w| — 



-Q=0 = -tS^T^ I dck I dx I dy I dqi / dq2 



UJm 1^771 j^'rrz 

1 1 



(iujm - e/c)(jw + iQm - - X - ql/2m) iuj'^ - eu + y - ql/2m 

1 1 



i^;, + ~ ek + x + y - {qi+ 172)^/ (2m) icj.;„ + - + x + y - {qi - q2)'^/2m 
Integrating over y first, and then over x, we obtain 



(F6) 



"Q=0 



7 (0) /■ /■ f 

__^j,3 y2 / / / dg2 [sgn(w„,. + f7,„) + sgnrj„i] [sgnw^ - sgn(tj^ + f7„)] 

7^ «fJ ,., o "'-00 Jo 

(F7) 



(16 

1 1 



it^m - Efc JW™ + 2ifi,„) - ek - qiiqi + q2)/m' 
Integrating next over ek, then over (72, and finally over qi using 



1 TTTTi 

rf<72— 7 ^ nr^^ ^TTT- = * ^g^^ [sgnwm - sgn(a;m + 2f]m)] , (F8) 

+ g2)/(TO) ~ 2zi2,m gi 



we obtain, to logarithmic accuracy 



-Q-o-~j^^TY:in^^Sn, (F9) 
where 

Sn^T^ [sgnwm - sgn(w™ + 20^)] [sgn(u;„ + 17^) + sgnO^] T ^ [sgnw^ - sgn(w^„ + fl^)] . (FIO) 
Performing the frequency summation, we finally obtain 

=^ - (Fll) 



(27r4) ™ 



where uq = mC/(0)/(27r). In Sec. lIVBI we have shown that 



^E"'1^W = -2^'C(3) + ..., (F12) 



where the dots stand for regular T^,T^, etc terms [see Eq. (|4.36fl and Eq. 14.37|l ]. Thus, 

2^3 



25=o = C(3)u^ — ■ (F13) 



The evaluation of the 2fcF-contribution from diagram 2(b) and of the entire contribution of diagram 2(a) (the latter 
involves only a process in which one of the vertices carries momentum near zero and the other one near 2k p) proceeds 
in the same way. The net result from the two diagrams is 

E"" = Ci3)iul + ul,^-uoU2k,)^. (F14) 

T^Vp 

In all three contributions, the total momentum in each vertex is near zero, i.e., both the vertices with momentum 
transfers near zero and near 2kp, have the structure of (k, — k; k, — k) and (k, — k;— k, k), respectively. Calculating 
the specific heat corresponding to S in Eq. (|F14|) . we find that the result for C{T) coincides with that in Eq. (|4.24() . 
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The remaining task is to show that the contribution from forward scattering, i.e., from processes of the type 
(k, k; k, k), vanishes. To see this, we repeat the same calculation, assuming now that the momenta in all four Green's 
functions in Fig. |3 2(b) are nearly parallel to each other. Instead of Eq. ljF6|l . we then obtain 



"Q=0 



2mL^ 3 



E 



dei- 



dx 

> 

1 



dy I dqi 



dq2 



{iLU - €k){i{ijJ + ri) - Efc - x - 9^/(2™)) 

1 



ek - X - y - {qi - q2)^/2m 



Integrating first over y, and then over x, we obtain 

1 



"Q=0 



d€k 



iujm - e/c \ii^m - Cfe + qiqi/m iujm - e/c - qiq2/m 



(F15) 



This integral vanishes as the poles of the integrand are located in the same half-plane of e^. This completes our proof 
of the statement that, to second order in the interaction, the singular part of S comes only from backscattering,i.e., 
from the diagrams containing vertices with zero total momentum. 



2. third order in the interaction 



Next, we consider the diagrams for ^ to third order in U . As we have shown in Sec. IIVBI the non-analyticity in 
results from the logarithmic singularity in the momentum integral, followed by the Matsubara sum 

rV / dQQ^ aTVr!f„ln|f7„| (xT^. 

At the second order, the singular term in the momentum integral, fJ^j/Q^ is obtained by multiplying the dynamic parts 
of two polarization bubbles. Higher-order diagrams contain higher powers of the polarization bubbles; however, we 
still need to select terms that behave only as /Q^: both less or more divergent terms do not result in a logarithmic 
non-analyticity of the momentum integral, and hence in a non-analyticity in S. At the third order, such terms are 
obtained by selecting the dynamic parts of two out of three polarization bubbles, while putting = in the third 
bubble. A subtle point here is that the static part of the bubble involves the integration over all internal fermion 
momenta, 1, i.e., in contrast to the dynamic part, there is no correlation between the directions of boson momentum, 
Q, and 1. Hence, vertices that appear in the third order diagrams are generally not "one-dimensional" in a sense 
that the four fermion momenta are not directed along the same line. Furthermore, despite the fact that the static 
polarization for free fermions in 2D n(ri = 0, Q) is a constant (= — m/27r) for all Q < 2kp, only H(il = 0, Q 0) comes 
from the states in the immediate vicinity of the Fermi surface. For generic Q, the static polarization bubble involves 
fermion states away from the Fermi surface [50l |. As a result, the convolution of U{Q) and two fermion propagators, 
which form a polarization bubble with non-zero external momentum, cannot be expressed as the angular average 
of bare interaction between the particles on the Fermi surface. As an example, consider the "Q = 0" contribution 
from diagram 3(b) in Fig.0 This contribution involves only bubbles with small external momenta, i.e., all internal 
fermion momenta-p, k, and 1-are located near the Fermi surface. Two of these momenta, e.g., p and k, must be 
nearly antiparallel; otherwise there is no non-analyticity in S. However, the direction of the third momentum-l-with 
respect to k is arbitrary. As a result, depending on the choice of two antiparallel momenta, the "Q = 0" contribution 
from diagram 3(b) contains a term proportional to UqUt^, where 

ug = {m/2n)Ui2kF sin6'/2), (F16) 

as well as another one proportional to u^{uq), where 

{ug) = - / u{e)de. (F17) 

In this last term, the integration goes over the entire Fermi surface. [It is understood that w^r = U2kp-] 

The same reasoning applies to the 2fci?-contribution from this diagram. In addition, the 2fci?-contribution involves 
the convolutions of the interaction potential with the Green's functions forming bubble. As we just said, the corre- 
sponding momentum integral is not confined to the Fermi surface. Evaluating the integrals, we find that the 2fci?-part 
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of diagram 3(b) contains a term proportional to u^uq and another one proportional to u^llug)), where 



1 r ./I + sin 61/2 + Jl-sin9/2 

{{ue)) ^- dOue cos - In ^ . . ^ . / (F18) 

^ -^/l + sm6'/2 — -^/l — smt'/z 



TT 



(These two terms occur for different choices of two anti-parallel momenta.) If u does not depend on 6, then {{ug)) — u. 
The first term -proportional to uJuq" involves the static bubble, 11(0, 2fci?), which, once again, is determined by the 
states far away from the Fermi surface. This term is similar to the "Q = 0" contribution from the states near the 
Fermi surface, simply because for free fcrmions in 13 = 2, Ii{2kp) ~ n(0). To distinguish between the Fermi-surface 
and non-Fermi surface contributions, we multiply 11(2^^) by ((1)) (according to Eq. I|F18I ((1)) = 1) to emphasize 
that the integration is not confined to the Fermi surface. 

Applying this reasoning to all third order diagram, and combining all choices of choosing two dynamic and one 
static bubble, we find 

Ssa = - (uo(w0U7r-e) + 2uou^{ue)) K + {u^((ug)) + 2uou^((u0))) K (F19a) 
S36 = {^ul{ue) + 2ulu^) K + [Aul{{ue)) + 2uluo{{l))] K (F19b) 
^^c = -^[ul + ul{{l))]K (F19c) 
= 2u^{ueu^-e)K + 2uo{{ul))K, (F19d) 

where K = C{3)T^ /nvj,. Adding up Eqs. HF19alF19dp . we obtain a total third-order contribution to S 

S3 [-4ul + 2ulu,, + {Aul- 2uou^){ue) ~ {uq -2u^){ugu^_g)] K 

+ [i-Aul + 2uoul){{l)) + {Aul - 2uou^){{ue)) - (u^ - 2uo){{uj))] K. (F20) 

Combining the last expression with the second-order result, Eq. I|F14|I . and using relation Eq. H4.2I) between S and 
C(T), we obtain for C{T)/T to third order in u 

SC (T) /T=- ^^^^ J- [(2^0 - + 2{{ul)) - (ueu^-e))^ + 3{u^ + {ueu^-e)f] , (F21) 

where 

UQ = UQ (1 -2uo + 2(ue)), ^ (1 - 2u^((l)) + 2((ue))) . (F22) 



Next, we verify whether Eq. (|F22|I can be obtained by substituting the renormalized static vertices into the second- 
order result for the specific heat, Eq. I|F14|I . To first order in u, we simply have: r'^(k, — k; k, — k) — uq and 
r'^(k, — k; — k, k) — U2kF- To evaluate the third-order contribution to C{T), we need to renormalize the vertices 
up to second order-the third order terms will then result as cross-products of first and second-order terms. Evaluat- 
ing the vertex corrections, presented diagrammatically in Fig. |31 in the same way as we evaluated the diagrams for 
S, we obtain 

r'=(k,-k;k,-k) = ^ [no + {{ul))] ; (F23a) 

r'=(k,-k;-k,k) = — [u, + (ugu^^g)] . (F23b) 
m 

Replacing uq and U2kF in Eq. ljF14p by renormalized vertices, Eq. ljF23a|) and Eq. (|F23b|) . correspondingly, we find 
after simple manipulations that it does indeed reproduce Eq. ljF21|l . This proves, to order u^, that the non- analytic 
term in the specific heat is expressed in terms of renormalized static vertices r'^'(k, — k; k, — k) and r* (k, — k; k, — k), 
i.e., in terms of r'^(7r). 

Using now the relations between r(k, — k; k, — k) and r(k, — k; k, — k) and the spin and charge components of r'''(7r) 



r(fc, -k; k, -k) = r^(7r) - r,^(^), r(fc, -k- -k, k) - -2v1{n), (F24) 

and restoring quasi-particle Z factors and m* /m, which come from self-energy insertions not considered above, we 
obtain Eq. H4.52|l . 

It is also instructive to re-express F'^ (tt) not in terms of the bare interaction potential, but in terms of another vertex- 
T^{6), which, we remind, is the static vertex renormalized by the states away from the Fermi surface. To this end, we 
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separate the Fermi-surface and non- Fermi-surface contributions to Eq. ljF23al lF23b|l . i.e., re- write Eq. (jF23al lF23b|) 
to order as 

r'=(k,-k;k,-k) = — [uo + {{ul))][l-2uo + 2{ug)]; (F25a) 
m 

r^-(k, -k; -k, k) - —u^[l~ 2u^{{l)) + 2{{ug))) (1 + ^"'""-'^ )]. (F25b) 
m 

(For simplicity, we neglected the Z-factor and effective mass renormalizations here and in what follows.) The first 
brackets in both formulas come from the states away from the Fermi surface, i. e, they give F'^ (tt) . The second brackets 
come from states near the Fermi surface and account for the difference between r*"' and F". Introducing spin and 
charge components of F'^ in the same way as in Eq. (|F24I) . i.e., as 

r'=(fc, _fc, k, -k) = - [f,{n) - fs{n)] , r\k, ~k; -k, k) = -2- ^ 
m m 

F'^(fc,-fc,fc,-fc) = -[7,(7r)-7,W], F'^(fc,-fc;~fc,fc) = -2- 7,W, (F26) 
m m 

we obtain 

7jfW = u^-\uA^-'2u^m + ^{ug)))^{{ul)) 

7s"W - -iM.(l-27i.((l»+2((ue))). (F27) 

Substituting Eq. ljF27p into Eq. ljF25bl) and (jF24l) . we obtain the relation between fa and 7" presented in the main 
the text [ Eq. gSl]- 

Finally, when the interaction is strongly peaked at Q = 0, so that Uq is much larger than ug for a generic 9, including 
9 = TT, only corrections to uo matter. These corrections come from the ring diagrams and can be summed up exactly. 
The full combinatoric factor for the Uq term from the ring diagram of order n is (— l)"2"~^(n — 1). Evaluating the 
sum over n, we reproduce Eq. H4.56|l . 
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